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ABSTRACT 


Risk  assessment  of  aging  aircraft  components  can  be  achieved  by  operational  de¬ 
rating  using  a  safety  factor  subjectively  selected  from  experience  and  heuristics.  This 
investigation  involves  synthesizing  currently  available,  maturing  computer-aided  methods 
into  a  format  of  objective  quantitative  risk  assessment.  The  methodology  is  applied  to 
quantify  the  effect  of  corrosion  on  P-3C  main  landing  gear  lower  drag  struts.  This  kind  of 
synthesis  is  appropriate  wherever  structural  operational  risk  is  a  concern.  The  P-3  has 
undergone  many  modifications  since  the  1950’s  and  the  lower  drag  struts  are  being 
scrapped  due  to  internal  surface  corrosion.  The  corrosion  process  is  random,  resulting  in 
pits  varied  spatially  and  in  severity.  These  corrosion  attributes  are  merged  into  a  one 
random  variable  probability  model.  The  casual  relation  of  the  corrosion  to  structural  load 
is  analyzed  by  finite  elements.  The  structural  configuration  model  input  is  provided  by 
computer-aided  drafting,  verified  by  physical  measurement.  The  effect  of  corrosion  on 
current  strut  population  reliability,  as  well  as  the  future,  is  computed.  The  conclusion  is 
that  even  under  severe  corrosion,  compressive  buckling  is  not  an  issue.  All  the  other 
failure  modes  (compressive  yielding,  tensile  yielding,  and  fracture  by  fatigue)  can  be 
assured  by  one  cold  temperature  proof  test. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

The  P-3C  “Orion”  is  a  four  engine,  low-wing  aircraft  designed  for  patrol  and 
antisubmarine  warfare  operated  by  the  U.S.  Navy  and  several  other  countries.  The 
aircraft  was  originally  designed  and  built  by  the  Lockheed  Aircraft  Corporation  in 
Burbank,  California  beginning  in  the  late  1950s.  It  has  undergone  several  upgrades  and 
modifications,  including  clearance  to  operate  at  heavyweight  loading  conditions  in  1966. 
Figure  1-1  shows  a  P-3C  in  flight  with  a  full  load  of  armament.  The  landing  gear  are  in  a 
tricycle  arrangement  with  two  main  gear  and  one  nose  gear.  All  three  landing  gear 
incorporate  hollow  cylinder  components  that  are  open  to  the  environment  and  have 
experienced  varying  degrees  of  internal  surface  corrosion  after  many  years  of  fleet 


Figure  I- 1.  P3-C  “Orion”  Aircraft. 
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operation.  This  investigation  specifically  addresses  the  reliability  of  the  main  landing 
gear  lower  drag  struts,  which  have  experienced  the  most  severe  surface  corrosion  and  are 
currently  being  replaced  during  aircraft  rework  at  higher  than  expected  rates.  Finite 
element  analysis  (FEA),  with  a  probabilistic  approach,  and  identification  of  a  direct 
filtering  technique  for  weak  struts  were  combined  to  provide  the  decision  maker  with  a 
quantitative  risk  assessment  tool. 

B.  LANDING  GEAR  DRAG  STRUT 

Each  P-3C  landing  gear  consists  of  dual  wheels  with  forward  retracting  struts,  as 
shown  in  Figure  1-2.  The  main  struts  on  each  of  the  three  gear  support  the  weight  of  the 
aircraft.  The  drag  struts  hold  the  main  struts  in  place  by  carrying  transverse  loads  from 
conditions  such  as  wheel  spin-up  during  landing  and  braking  in  both  forward  and  reverse 
directions.  The  drag  struts  are  hollow,  tapered  cylinders  made  from  AISI 4340  steel 
forgings,  flash  welded  together  in  the  middle,  as  shown  in  Figure  1-3.  There  are  small 
holes  left  in  the  ends  of  the  enclosed  portion  of  the  strut  so  that  hot  gases  can  escape 
during  the  flash  welding  process.  The  tapered  ends  are  machined  into  lugs  that  hold  balls 
for  attachment  to  neighboring  hardware  thus  putting  the  drag  strut  in  a  pin-ended 
boundary  condition. 
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Figure  1-3.  Main  Landing  Gear  Lower  Drag  Strut. 


The  P-3C  is  currently  being  overhauled  under  the  Sustained  Readiness  Program 
(SRP)  to  improve  the  reliability  of  the  aircraft  to  allow  it  to  safely  fly  up  to  the  specified 
original  service  life  period.  The  main  landing  gear  lower  drag  strut  is  one  of  many 
components  that  are  inspected  during  SRP.  One  of  the  drag  struts  that  was  rejected 
during  this  inspection  process  was  sent  to  the  Naval  Postgraduate  School  (NPS)  for  use  in 
this  investigation.  This  strut,  referred  to  as  the  sample  strut  in  the  remainder  of  this 
report,  was  evaluated  at  NPS  by  the  author  to  extract  required  data  using  the  detailed 
procedures  described  in  Appendix  A. 

C.  LOADING 


The  various  takeoff,  landing,  and  taxi  operations  that  the  P-3C  experiences  in 
operation  place  the  drag  strut  in  both  tension  and  compression.  The  maximum  tension 
condition  occurs  during  two  wheeled  braking  at  high  speed  with  a  defined  limit  load  of 
160,571  lbs.  The  maximum  compression  condition  occurs  during  two  wheeled  braking 


while  the  aircraft  is  backing  down  following  the  use  of  reverse  thrust,  with  a  defined  limit 
load  of  132,589  lbs.  The  likelihood  of  occurrence  of  this  condition  is  extremely  remote 
because  the  Navy  trains  its  pilots  not  to  use  wheel  brakes  to  stop  the  aircraft  while 
backing  down.  The  next  highest  compression  load  occurs  during  landing  when  the 
landing  gear  spring  back  following  initial  contact  with  the  runway.  This  limit  load  is 
defined  as  80,033  lbs.  These  loading  conditions  are  for  a  heavyweight  aircraft  and  were 
obtained  from  Lockheed  Report  (LR)  13680.  [Ref.  1] 

D.  FAILURE  MODES 

There  are  three  basic  ways,  or  modes,  in  which  the  drag  strut  could  potentially  fail 
in  service.  (1)  Under  the  tension  loading  condition,  the  strut  could  fracture  due  to  fatigue 
cracking  or  manufacturing  defect.  This  would  be  a  catastrophic  failure  that  would  cause 
collapse  of  the  main  landing  gear.  (2)  Under  the  compression  loading  condition  the  strut 
could  experience  column  buckling,  followed  by  gross  plastic  deformation.  This  sort  of 
damage  would  more  than  likely  cause  only  partial  collapse  of  the  landing  gear  because 
there  would  still  be  a  physical  connection  between  the  components.  (3)  The  strut  could 
simply  plastically  yield  in  a  way  that  the  landing  gear  geometry  would  be  altered  such 
that  the  retraction,  extension,  or  locking  processes  would  not  function  properly.  This 
investigation  specifically  addresses  the  column  buckling  failure  mode  by  analysis  and 
suggests  direct  filtering  by  low  temperature  proof  test  to  address  the  fracture  and  plastic 
yielding  modes. 

E.  INTERNAL  CORROSION 

Because  of  the  pressure  relief  holes  in  the  ends  of  the  strut,  the  interior  is  open  to 
the  environment  and  subjected  to  atmospheric  conditions  as  the  aircraft  climbs  and 
descends.  Moist  air  condenses  on  the  interior  surfaces  of  the  strut  and  eventually  causes 
surface  corrosion.  The  original  Lockheed  design  included  a  synthetic  rubber  coating  on 
the  interior  surface  to  prevent  this  sort  of  corrosion.  Over  the  years,  several  depot  or 
factory  level  preventative  maintenance  procedures  have  been  incorporated  to  solve  other 
structural  reliability  issues,  including  hydrogen  embrittlement  and  surface  stress  relief 
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following  machining.  These  procedures  call  for  baking  the  strut  at  temperatures  of  up  to 
400  degrees  Fahrenheit  which  destroys  the  synthetic  rubber  coating  on  the  interior  [Ref. 
2].  Once  the  coating  is  destroyed,  the  accumulation  of  moisture  provides  an  environment 
conducive  to  corrosion. 

F.  CORROSION  CHARACTERIZATION 

The  corrosion  found  in  the  sample  rejected  drag  stmt  had  several  notable 
characteristics.  An  example  of  a  semi-cylindrical  section  of  the  stmt  with  flaking  surface 
corrosion  is  shown  in  Figure  1-4.  Once  the  scale  was  removed,  the  pitting  found  on  the 
surface  was  generally  grouped  in  three  spatial  locations  on  the  thin  walled,  cylindrical, 
middle  section  of  the  stmt.  The  three  locations,  indicated  in  Figure  1-5,  corresponded  to 
regions  where  moisture  could  naturally  accumulate  during  different  phases  of  operation 
on  the  aircraft.  These  regions  would  either  be  on  the  top  or  bottom  of  the  cylindrical 
section  depending  on  whether  the  landing  gear  were  extended  or  retracted.  The  observed 
pattern  was  consistent  with  the  general  understanding  of  the  corrosion  process  as 
described  in  Appendix  B.  Specifically,  the  observed  locations  were  regions  where  water 
would  pool  on  the  lower  surface  or  water  vapor  would  condense  on  the  upper  surface.  In 
the  corrosion  process,  the  presence  and  longevity  of  moisture  on  a  surface  is  directly 
related  to  the  level  of  subsequent  corrosion.  There  was  virtually  no  corrosion  on  the 
regions  that  would  be  on  the  sides  of  the  cylindrical  section. 

The  small  pits  were  generally  semi-spherical  in  shape  to  the  naked  eye  with 
highly  irregular  boundaries  when  viewed  through  a  microscope.  The  larger  pits  appeared 
to  have  been  formed  when  the  smaller  spherical  pits  grew  together.  Their  irregular 
boundaries  were  visible  without  the  aid  of  a  microscope  and  appeared  to  be  influenced  by 
the  molecular  imperfections  in  the  material  itself.  Some  of  the  pits  were  broad  and 
shallow  while  others  where  very  narrow  and  deep.  Quantifiable  characteristics  of  the 
corrosion  including,  pit  depth,  pit  volume,  and  pit  location  were  measured  and  are 
included  in  Appendix  A. 
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Figure  1-4.  Semi-Cylindrical  Section  With  Flaking  Corrosion. 
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Figure  1-5.  Regions  of  Corrosive  Pitting  Concentration. 


G.  CURRENT  INSPECTION  TECHNIQUES 

At  the  time  of  this  investigation,  the  Navy  was  not  performing  an  active 
inspection  of  fleet  aircraft  to  determine  the  condition  of  drag  struts  in  service.  An 
inspection  technique  had  been  developed  that  included  X-ray  to  identify  areas  of 
corrosion  followed  by  ultrasonic  inspection  to  determine  wall  thickness.  Even  by  the  best 
operators  with  the  best  equipment  these  methods  are  indirect,  in  that  they  do  not  measure 
the  strength  of  the  strut.  They  characterize  the  effect  of  corrosion  on  the  strut,  to  a 
limited  extent.  A  criterion  had  been  established,  based  upon  a  recommendation  from  the 
manufacturer,  that  any  strut  with  pitting  that  resulted  in  a  wall  thickness  less  than  0.140 
inches  should  be  rejected  [Ref.  3].  Examination  of  the  X-rays  that  were  taken  of  the 
sample  strut,  and  information  received  from  the  Navy’s  representative  at  the  SRP 
contractor’s  plant,  showed  that  these  inspection  techniques  were  highly  unreliable  for  this 
situation.  The  author  was  told  that  the  X-ray  and  ultrasonic  inspections  on  the  sample 
strut  did  not  reveal  pitting  in  excess  of  the  current  rejection  criteria;  however,  the  strut 
was  rejected  on  a  “gut  feel”  because  a  “very  unusual  pit”  was  detected.  Dissection  and 
examination  of  the  sample  strut  revealed  that  at  least  34  pits  existed  that  were  in  excess  of 
the  current  rejection  criteria.  The  minimum  wall  thickness  was  found  to  be  0.122  inches. 

H.  ANALYSIS  METHODOLOGY 

The  methodology  used  in  this  analysis  has  been  summarized  in  Figure  1-6  which 
is  a  flow  chart  of  the  steps  necessary  to  arrive  at  the  ultimate  risk  assessment. 
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II.  PROBABILISTIC  RELIABILITY 


A.  CONTINUOUS  RANDOM  VARIABLES 

To  determine  the  reliability  of  a  component  such  as  the  drag  strut  one  must 
recognize  that  there  is  a  finite  population  of  drag  struts  in  existence.  The  population  has 
distinct  characteristics,  or  properties,  in  terms  of  many  different  variables  such  as 
strength,  weight,  physical  dimensions,  levels  of  corrosion,  etcetera.  These  variables  are 
random,  or  stochastic,  because  there  is  no  predetermined  order  or  exact  set  of  values  for 
the  variables  in  any  given  strut.  The  variables  are  continuous  because  they  can  take  on 
any  value  within  certain  bounds  described  by  the  geometry,  materials,  and  manufacturing 
processes  that  make  up  the  strut.  Each  individual  drag  strut  has  its  own  levels  of  these 
continuous  random  variables  (RV).  There  are  mean,  or  average,  values  for  an  individual 
strut  but  different  values  at  any  given  location  within  the  strut.  There  is  no  way  to  predict 
the  exact  levels  of  these  variables;  however,  one  can  put  bounds  on  the  levels  given  some 
knowledge  about  the  population.  More  importantly,  one  can  establish  the  probability  of 
finding  a  particular  strut  with  a  given  level  of  one  of  these  variables  if  the  probability 
distribution  function  for  the  population  is  known. 

B.  PROBABILITY  DISTRIBUTIONS 

Let  the  probability  of  finding  a  discrete  value  x,  of  some  random  variable  X,  be 
given  by  P(X  =x)=  f(x),  where  f(x)  is  called  the  probability  distribution  function.  In 
probability  distributions  of  continuous  RVs,  the  probability  that  a  given  variable  X  takes 
on  any  one  discrete  value  is  generally  zero.  However,  one  can  determine  the  probability 
that  X  lies  between  two  different  values  and  thus  the  probability  distribution  for  a 
continuous  RV  can  be  found.  A  real-valued  RV  must  definitely  lie  between  and  °o, 
and  the  sum  of  all  the  probabilities  must  equal  unity,  because  the  entire  population  is 
being  considered.  These  statements  are  represented  by: 

f(x)>  0  (1) 
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(2) 


f*  f(x)dx  =  1 

J—O O 

P(a<X  <b)=[b  f(x)dx  (3) 

The  probability  that  the  random  variable  X  takes  on  at  least  some  value  jc  is  given 
by  the  cumulative  distribution  function  (CDF)  F(x),  where 

F(x)=  P(X  <  x)=  P(-oc<X<  x)=  j\f(u)du  (4) 

C.  DISCRETE  RANDOM  VARIABLES 

Continuous  RVs  are  fine  for  mathematical  relationships;  however  one  can  only 
measure  or  quantify  discrete  values  of  random  variables  when  dealing  with  a  real  object. 
The  discrete  values  are  dependent  upon  the  scale,  or  resolution,  of  the  instrument  used  to 
measure  or  evaluate  the  object.  If  the  scale  is  small  enough,  then  the  attributes  of  the 
discrete  RV  approach  that  of  the  continuous.  In  this  case,  the  summation  of  all  of  the 
measured  discrete  values  of  the  RV  must  be  unity,  because  the  entire  population  is  being 
considered.  This  is  represented  by: 

£/W=  1  (5) 

j: 

The  remainder  of  this  investigation  deals  with  discrete  RVs.  In  each  case,  the 
scale  was  chosen  as  a  compromise  between  physical  limitations  of  measuring  tools,  the 
time  available  to  quantify  the  RV,  and  a  desire  to  represent  the  real  continuous  RV  as 
closely  as  possible. 

D.  SPECIFIC  DISCRETIZED  RANDOM  VARIABLES 

The  drag  strut,  as  originally  manufactured,  was  thought  to  be  reliable  enough  for 
the  planned  service  life  of  the  aircraft.  To  date,  there  have  reportedly  been  no  in-service 
structural  failures,  which  validates  the  original  analysis.  However,  that  analysis  did  not 
consider  the  effect  of  corrosion  because  that  damage  process  had  been  precluded  by  the 
synthetic  rubber  coating.  The  RVs  that  were  important  for  determining  the  original  strut 
reliability  were  controlled  during  the  manufacturing  process.  The  material  had  to  meet 
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certain  specifications  in  terms  of  hardness,  modulus  of  elasticity,  ultimate  strength, 
etcetera.  The  stmt  had  to  be  machined  to  within  certain  dimensional  tolerances  and  the 
flash  weld  had  to  be  performed  under  specific  process  specifications.  Finally  the 
manufactured  stmt  was  proof  tested  in  tension  at  room  temperature,  which  served  as  a 
high  pass  filtering  process  to  remove  any  weak  stmts  that  slipped  by  the  other  control 
mechanisms.  The  proof  test  load  was  set  at  67  percent  of  the  material  yield  strength  or 
220,600  pounds  and  was  specifically  intended  to  assure  the  reliability  of  the  flash  weld 
[Ref.  4].  The  RVs  that  are  important  for  determining  corroded  stmt  reliability  have  now 
been  augmented  by  things  relevant  to  the  corrosion  damage  process.  They  differ  from 
the  original  RVs  in  that  their  levels  can  not  be  controlled,  but  only  evaluated.  The 
augmented  RVs,  or  corrosion  induced  effects,  can  be  characterized  by  either  the  spatial 
location  of  one  or  more  damage  sites  on  the  stmt,  or  the  individual  attributes  of  any  one 
damage  site.  These  damage  sites  will  be  referred  to  as  pits  from  this  point  forward. 

1.  Pit  Spatial  Location 

a.  Single  Pit 

There  are  infinitely  many  spatial  locations  for  pits  within  the  stmt  because 
location  is  a  continuous  RV.  Location  had  to  be  converted  to  a  discrete  RV  for 
characterization  and  analysis.  The  critical  part  of  the  stmt  for  resistance  to  column 
buckling  is  the  thin  walled  cylindrical  portion  in  the  middle.  This  portion  was  divided 
into  four  semi-cylindrical  sections  as  shown  in  Appendix  A.  These  sections  were  further 
subdivided  into  an  18  by  24  grid  that  contained  432  possible  locations.  These  locations 
will  be  referred  to  as  grid  elements  or  just  grids  in  the  remainder  of  this  document.  But 
column  buckling  is  particularly  sensitive  to  imperfections  at  critical  locations; 
specifically  the  anti-nodes  of  the  harmonic  modes.  For  this  reason,  special  emphasis  was 
given  to  cases  where  the  worst  pits  were  located  at  critical  nodes.  This  reduced  the 
number  of  combinations  to  be  analyzed  to  a  manageable  level. 
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b.  Pit  Clustering  or  Multiple  Pits 

When  considering  multiple  pits  within  the  discrete  location  RV,  there  are 
still  more  possible  combinations  of  pits  than  could  be  evaluated  in  this  investigation. 
Because  the  order  in  which  the  locations  were  selected  was  not  important,  the  number  of 
combinations  possible  could  be  calculated  from  the  following  formula  for  combinations 
of  n  things  taken  r  at  a  time: 

=  'v;  :  n\=n(n-l)(n-2}~l  (6) 

r\{n-r)l 

If  two  pits  were  considered,  there  would  be  93,096  combinations  of  pit 
spatial  location.  For  three  pits,  there  would  be  over  13  million  combinations.  These 
numbers  are  for  only  one  kind  of  pit.  If  pits  of  different  characteristics  were  coupled 
together  with  the  possible  location  combinations,  the  numbers  would  get  even  larger.  To 
deal  with  this,  a  series  of  extreme  cases  were  chosen.  The  analysis  progressed  from  pits 
to  through-the-thickness  holes  of  increasing  diameter,  up  to  a  maximum  of  two  inches. 
The  holes  were  assumed  to  be  a  practical  determination  of  the  worst  possible  case,  i.e.,  no 
operator  would  miss  or  accept  a  two  inch  hole  in  the  side  of  any  strut. 

2.  Pit  Attributes 

There  are  multiple  attributes  that  could  be  used  to  describe  the  corrosive  pits 
found  within  the  drag  strut.  Each  of  them  is  a  continuous  RV  in  its  own  right  which  has 
to  be  discretized  in  order  to  be  used  in  analysis.  The  following  attributes  were  chosen  on 
the  basis  that  they  would  have  a  distinct  impact  on  the  strength  or  reliability  of  the  strut. 
The  detailed  procedures  used  to  extract  the  values  of  these  parameters  in  the  sample  drag 
strut  can  be  found  in  Appendix  A. 

a.  Pit  Depth 

This  was  the  easiest  attribute  to  quantify  because  it  could  be  measured 
directly.  There  was  no  special  technique  required  to  evaluate  the  impact  of  pit  depth.  In 
other  words,  a  deep  pit  is  worse  than  a  shallow  pit  because  there  is  a  direct  reduction  in 
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the  amount  of  material  present  to  carry  stress  within  that  grid.  Each  grid  element  was 
discretized  by  the  depth  of  the  deepest  pit  within  that  grid. 

b.  Pit  Volume 

The  volume  of  material  missing  due  to  corrosion  was  more  difficult  to 
quantify  because  it  could  not  be  measured  directly.  Instead,  the  amount  of  liquid  required 
to  fill  the  pits  in  each  grid  was  measured.  Alcohol  was  chosen  for  this  because  of  its  low 
surface  tension  which  allowed  it  to  spread  evenly  across  the  pits  and  not  bead  up.  This 
attribute  should  also  directly  relate  to  the  strength  of  the  strut  because  of  the  reduction  in 
material  present  to  carry  stress  and  provide  stability. 

c.  Pit  Shape 

The  shape  of  individual  pits  would  have  a  significant  impact  on  a  fracture 
mechanics  or  fatigue  life  analysis.  However,  column  buckling  is  more  sensitive  to  the 
overall  cross-section  moment  of  inertia  and  the  overall  distribution  of  material,  not  the 
shape  of  a  single  pit  on  the  interior  surface.  Because  the  fracture  and  fatigue  failure 
modes  are  being  addressed  by  proof  test,  except  for  observing  the  general  nature  of  pit 
shape  and  growth  as  described  above,  there  was  no  attempt  to  record  the  pit  shape  in  each 
grid. 


d.  Equivalent  Spherical  Pit  Radius 

One  goal  of  the  pit  attribute  definition  process  was  to  provide  a  means  for 
modeling  corrosion  with  finite  elements  (FE).  To  this  end,  a  single  parameter  was 
required  to  combine  both  the  pit  depth  and  pit  volume.  The  chosen  parameter  was 
equivalent  spherical  pit  radius.  This  was  the  radius  of  a  sphere  such  that  a  slice  from  that 
sphere,  at  a  depth  equal  to  the  grid  pit  depth,  would  occupy  the  same  volume  as  the  grid 
pit  volume.  This  is  more  clearly  illustrated  in  Figure  II- 1,  where  R  is  the  equivalent 
spherical  pit  radius,  h  is  the  pit  depth  and  V  is  the  pit  volume.  The  justification  for  the 
choice  of  this  attribute  in  terms  of  corrosion  pit  propagation  is  included  in  Appendix  B. 
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Spherical  Radius,  R 


Pit  Depth,  h^j  ]  ^Pit  Volume,  V 

Cylindrical  Wall 
Figure  II- 1.  Equivalent  Spherical  Pit  Radius. 

e.  Pit  Severity 

The  final  attribute  used  to  quantify  corrosion  in  the  strut  was  defined  as  pit 
severity,  which  is  simply  the  grid  pit  depth  multiplied  by  the  grid  pit  volume.  It  is 
somewhat  like  using  the  first  moment  of  mass  to  find  the  centroid  of  a  solid.  The 
contribution  of  an  individual  element  to  the  centroid  of  the  object  depends  upon  its  mass 
and  its  distance  from  the  origin.  In  this  case,  the  contribution  of  pitting  in  an  individual 
grid  to  the  strength  of  the  strut  depends  upon  the  volume  of  material  missing  and  its 
distance  from  the  surface.  This  attribute  was  needed  so  that  a  probability  distribution 
could  be  developed  that  would  accurately  describe  the  chances  of  finding  a  certain  level 
of  corrosion  in  a  certain  location.  The  assumption  was  that  the  statistics  of  finding 
corrosion  at  any  grid  location  (for  the  sample  strut)  could  be  generalized  to  be  the 
probability  of  occurrence  of  corrosion  severity  at  any  location,  within  any  strut,  from  any 
P-3  aircraft.  This  thought  is  based  upon  fractal  theory;  or  the  idea  that  within  a 
population  of  drag  struts,  the  map  of  the  corrosion  severity  within  a  strut  can  be 
duplicated  in  any  other  strut  from  the  population.  All  that  is  required  is  the  appropriate 
scale  for  the  map,  or  in  this  case  the  grid.  The  justification  is  that  the  sample  strut  came 
from  the  population  subset  of  rejected  struts.  It  contains  the  character  of  the  worst 
(because  it  was  rejected)  and  there  is  a  distinct  distribution  of  corrosion  (not  uniform 
corrosion  all  over). 
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E.  PIT  SEVERITY  PROBABILITY  DISTRIBUTION 

The  pit  severity  at  each  grid  location  was  plotted  as  a  3D  surface  as  shown  in 
Figure  II-2.  The  plot  has  been  oriented  to  allow  viewing  orthogonal  to  the  xy  plane 
(location  axes),  with  10  levels  of  color  providing  discrimination  along  the  z  axis  (pit 
severity).  Each  of  the  432  grid  boxes  was  further  divided  into  fifths  to  provide  additional 
resolution  in  location.  The  number  of  fifth  grids  at  each  of  the  10  levels  of  pit  severity 
was  counted  and  then  normalized  by  the  total  number  of  fifth  grids  or  2,160.  The  result 
is  the  frequency  of  occurrence  of  all  the  different  levels  of  severity  within  the  cylindrical 
portion  of  the  strut.  Because  this  frequency  of  occurrence  has  been  normalized  by  the 
total  number  of  fifth  grid  elements,  it  also  meets  the  requirements  of  a  probability 
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Figure  II-2.  Pit  Severity  Distribution  Surface  Plot. 
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distribution.  The  result  is  a  probability  distribution  of  corrosion  severity  in  the 
cylindrical  section  of  the  strut.  Recognizing  that  a  multiplicative  relationship  exists  in  the 
definition  of  pit  severity  and  the  corrosion  process  in  general  the  frequency  of  occurrence 
data  was  plotted  in  natural  log  space.  This  allowed  determination  of  the  mean  and 
standard  deviation  for  a  normal  distribution,  which  could  then  be  transformed  back  to 
regular  severity  space.  The  result  would  be  a  lognormal  distribution  of  the  probability  of 
finding  a  pit  of  a  certain  severity  in  any  given  drag  strut.  Figures  II-3  and  II-4  show  the 
results  of  this  process,  including  the  equations  for  the  respective  distributions.  The  ten 
individual  levels  of  severity  are  plotted  as  circles,  representing  a  histogram  of  all  the  grid 
severity  measurements.  The  reason  there  are  no  points  on  the  left  of  the  normal 
distribution  in  the  natural  log  space  is  that  pits  with  severity  at  that  level  would  be  at  the 
microscopic  level,  as  indicated  by  the  relative  size  icons  at  the  bottom  of  Figure  II-3. 


Figure  II-3.  Fractional  Occurrence  of  Severity  in  Natural  Log  Space. 
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Figure  II-4.  Fractional  Occurrence  of  Severity. 


The  significance  of  this  distribution  is  in  the  probability  of  finding  the  worst  pit  in 
any  strut.  Over  time,  the  right  tail  of  the  distribution  will  shift  to  the  right,  i.e.  the  rare  or 
worst  pits  will  get  worse.  The  area  on  the  left  (small  or  common  pits)  will  be 
redistributed  as  the  overall  level  of  corrosion  increases,  but  the  probability  of  occurrence 
of  the  rare  pit  will  basically  stay  the  same.  Based  on  the  fractional  measurement  of  the 
current  sample  strut,  the  probability  of  finding  the  worst  pit,  or  P(S),  was  determined  to 
be  0.00093.  In  addition  to  the  lognormal  distribution,  a  power  fit  relationship  was 
determined  as  an  estimate  of  the  upper  tail  probabilities,  because  the  lognormal  clearly 
does  not  sufficiently  model  what  was  observed.  Because  a  power  fit  has  been  used,  a 
more  clear  representation  of  the  function,  and  the  required  extrapolation  to  far  greater 
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levels  of  severity  can  be  seen  using  a  log-log  scale,  as  is  shown  in  Figure  II-5.  The  power 
fit  has  been  extrapolated  and  could  now  be  used  to  determine  the  probability  of  finding 
the  extreme  case  (worse  corrosion,  but  less  frequency  of  occurrence). 


Figure  II-5.  Extrapolation  of  Fractional  Occurrence  of  Severity. 

The  next  step  was  to  determine  the  probability  of  finding  a  pit  at  any  given 
location  within  the  strut.  The  pitting  was  assumed  to  be  independent  and  identically 
distributed  (iid),  i.e.,  the  probability  of  occurrence  of  a  pit  at  a  certain  location  in  the  strut 
did  not  depend  on  any  of  the  other  pits  in  the  strut  and  pitting  was  equally  likely  to  occur 
at  any  location.  This  actually  added  some  conservatism  to  the  analysis  because  the 
pitting  distribution  in  the  sample  strut  was  not  observed  to  be  equally  likely;  specifically 
there  was  no  pitting  observed  in  critical  locations.  However  this  assumption  was 
necessary  to  enable  the  continuation  of  the  probability  model.  Because  there  were  432 
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equally  likely  grid  elements,  then  the  probability  of  finding  a  pit  in  any  grid,  or  P(L),  was 
or  0.0023. 

Finally,  the  probability  of  finding  both  the  worst  pit  and  of  finding  it  at  a  specific 
location  was  the  joint  probability  or  P{S  C\L).  A  joint  probability  can  be  represented  by: 

P(A  nB)=  P(A)  P(B  I  A)  (7) 

where  P(B\A)  is  the  conditional  probability,  or  the  probability  of  B  occurring  given  that  A 
has  already  occurred.  When  the  “iid”  assumption  is  made,  the  conditional  probability 
simply  becomes  the  P(B)  and  the  joint  probability  becomes  the  product  of  the  two.  Thus 
P(S  n  L)  becomes  the  product  of  P(S)  and  P(L)  or  2.15  x  10"6. 

The  crucial  question  in  determining  the  column  buckling  reliability  of  the  drag 
strut  in  the  presence  of  corrosion  was:  What  is  the  probability  of  finding  a  pit  bad 
enough,  and  in  a  critical  location,  to  cause  buckling  in  any  given  strut?  The  answer  to 
this  question  was  formulated  through  realization  of  a  single  random  variable  by  finite 
elements.  Models  of  far  more  severe  corrosion  than  observed  in  the  sample  strut  were 
created,  and  the  probability  of  occurrence  of  those  configurations  were  determined  from 
the  joint  probability  using  the  product  of  the  power  fit  of  probability  of  severity  and  the 
probability  of  location  as  discussed  above.  This  physically  enormously  complex 
multivariate  problem  has  been  reduced  to  a  single  characteristic  random  variable 
encapsulating  the  effect  on  buckling  by  the  joint  attributes  of  all  the  variables. 

F.  RANDOM  VARIABLE  REALIZATION  BY  FINITE  ELEMENTS 

The  key  to  doing  a  risk  assessment  of  the  strut  in  buckling  was  to  be  able  to  relate 
the  analysis  to  a  random  variable  that  had  meaning  to  the  operator  or  decision  maker. 

The  envelope  within  which  the  aircraft  is  operated  is  defined  by  various  limitations. 
These  limitations  may  be  in  terms  of  airspeeds,  accelerations,  landing  sink  rates,  etcetera. 
The  manufacturer  had  done  the  work  of  translating  these  limitations  into  service  loads 
(forces  or  moments)  on  individual  components,  in  this  case  the  drag  strut.  Therefore,  all 
that  was  required  was  a  realized  random  variable  with  a  probability  distribution  that  was 
in  terms  of  load,  which  in  turn  can  be  related  back  to  operational  conditions. 
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There  were  at  least  two  ways  to  accumulate  the  statistics  of  a  realized  random 
variable  in  terms  of  load:  (1)  by  actual  buckling  testing  of  real  struts  or  (2)  by  finite 
element  analysis.  In  this  investigation,  FEA  was  used  to  perform  a  virtual  realization  of  a 
random  variable  in  terms  of  critical  buckling  load.  This  was  done  by  mapping  all  of  the 
attributes  described  above  into  a  single  random  variable  (which  is  a  measurable  quantity) 
and  then  finding  the  strength  (the  realized  random  variable)  by  finite  elements.  This 
made  it  possible  to  relate  the  service  load  to  the  critical  buckling  load  in  probability  terms 
and  hence,  perform  a  risk  assessment  as  will  be  presented  later  in  this  report. 
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III.  APPLICATION  OF  FINITE  ELEMENTS 


A.  ASSURANCE  OF  COMPUTATIONAL  RESULTS 

In  the  beginning  of  this  investigation,  certain  assumptions  were  made  concerning 
the  maturity  of  finite  element  analysis.  The  goal  was  to  provide  the  decision  maker  some 
analysis  results  that  could  be  taken  at  face  value,  with  an  obvious  visualization  of  the 
problem  and  its  solution.  Therefore,  the  investigation  was  guided  away  from  the  classical 
technique  reducing  the  problem  to  something  that  could  be  manually  calculated  through 
idealizations  and  assumptions.  In  that  case,  the  problem  is  pragmatically  idealized  to 
match  available  solution  techniques.  As  a  result,  the  solution  is  accurate  but  the 
interpretation  of  the  results  requires  a  great  deal  of  study  and  understanding  to  extend  the 
results  back  to  the  real  world.  Instead,  in  this  investigation  the  opposite  route  was  taken, 
that  is:  a  complex  model,  that  closely  matched  the  real  strut  was  constructed  and  a 
sophisticated  analysis  was  relied  upon  to  provide  results  that  could  be  directly  applied  to 
the  real  world. 

MSC/NASTRAN1  was  selected  because  it  was  a  well  tested  and  established 
software  application,  and  it  could  perform  the  automatic  solid  meshing  required  to 
represent  the  drag  strut  with  finite  elements.  NASTRAN  was  used  as  an  application 
package;  that  is,  no  effort  was  expended  in  element  development.  A  more  simple 
geometric  model  of  the  cylindrical  portion  of  the  strut  could  have  been  used,  but  the 
analysis  would  have  required  additional  compromise  in  forming  the  equivalent  boundary 
conditions.  The  interpretation  of  results  would  be  ambiguous  because  the  simpler  model 
would  need  to  be  related  back  to  the  actual  strut.  This  approach  was  not  used.  Instead, 
extensive  effort  was  expended  to  assure  the  geometric  model  was  a  duplication  of  reality; 
thereby  enabling  the  use  of  actual  boundary  conditions  and  results  that  could  be  used 
without  interpretation.  The  high  degree  of  replication  of  the  physical  drag  strut  by 
geometric  model  can  be  seen  in  the  modeling  section  to  follow,  say  Figure  IV- 1.  With 

1  MSC/  is  a  registered  trademark  and  service  mark  of  The  MacNeal-Schwendler  Corporation. 
NASTRAN  is  a  registered  trademark  of  the  National  Aeronautics  and  Space  Administration. 
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that  in  mind,  the  following  paragraphs  explore  some  of  the  details  of  FEA  that  are 
germane  to  the  results  obtained  and  provide  a  level  of  assurance  that  the  results  are  valid. 

B.  LINEAR  ELASTIC  BUCKLING 

The  analysis  of  the  drag  strut  in  this  investigation  was  performed  assuming  linear 
elastic  behavior  of  the  material  under  load.  In  linear  static  analysis,  a  state  of  stable 
equilibrium  is  assumed  at  all  times.  The  structure  deforms  under  load  in  a  linear  fashion, 
and  the  displacements  are  determined  directly  from  the  properties  of  the  material  being 
modeled;  specifically  the  elastic  modulus.  When  the  load  is  removed,  the  structure  is 
assumed  to  return  to  the  undeformed  position.  In  the  case  of  linear  buckling  analysis,  the 
solution  process  continues  until  the  assumption  of  stable  equilibrium  is  no  longer 
satisfied.  In  other  words,  at  some  point  the  loading  condition  will  be  such  that  the 
deformation  increases  without  any  increase  in  load.  The  structure  has  become  unstable 
and  has  reached  the  bifurcation  point.  It  is  assumed  that  there  is  no  yielding  of  the 
structure  up  to  this  point  and  the  direction  of  the  applied  load  has  not  changed. 

The  key  to  finite  element  analysis  is  the  stiffness  matrix.  The  stiffness  matrix  is 
determined  from  the  geometry  of  the  model  and  the  properties  of  the  elements  that  make 
up  the  model.  Included  in  the  properties  of  the  elements  are  the  characteristics  of  the 
material  that  the  element  is  representing,  in  this  case  the  modulus  of  elasticity  and 
Poisson  ratio.  The  size  of  the  stiffness  matrix  is  determined  by  the  complexity  of  the 
model  (the  number  of  nodes)  and  the  degrees  of  freedom  allowed  for  each  node.  The 
stiffness  matrix  is  the  minimum  potential  energy  of  all  the  nodal  displacements  in  static 
equilibrium.  Elastic  buckling  analysis  in  NASTRAN  includes  the  effect  of  differential 
stiffness.  This  allows  for  higher  order  strain  displacement  relationships  that  are  functions 
of  geometry,  element  type,  and  applied  loads.  This  can  be  viewed  as  a  linear 
approximation  of  softening  the  stiffness  matrix  under  compressive  load.  It  is  this 
softening  effect  that  leads  to  the  point  of  instability  in  the  structure. 
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C.  EIGENVALUE  DECOMPOSITION 

The  finite  element  solution  for  linear  buckling  is  determined  by  finding  the 
eigenvalues  of  the  stiffness  matrix.  The  eigenvalues  are  scale  factors  that  multiply  the 
applied  load  in  order  to  produce  the  critical  buckling  load.  The  eigenvalue  is  the 
characteristic  value  of  the  stiffness  matrix  such  that  the  effect  of  the  differential  stiffness 
matrix  has  caused  the  stability  of  the  structure  to  reach  zero.  This  is  the  point  of 
instability.  Mathematically,  the  eigenvalue  is  the  value  that  causes  the  determinant  of  the 
stiffness  matrix  to  go  to  zero.  The  process  of  finding  this  eigenvalue  is  called 
decomposition  in  matrix  algebra  terms.  There  are  many  eigenvalues  that  will  result  in 
this  point  of  instability;  as  many  as  the  number  of  rows  or  columns  in  the  stiffness  matrix 
since  it  is  a  square  matrix.  Generally,  only  the  smallest  eigenvalue  is  of  interest  because 
this  is  the  lowest  load  at  which  the  structure  will  become  unstable.  The  buckling 
eigenvalue  problem  reduces  to: 

[K  +  A,Kd]=0  (g) 

where  K  is  the  stiffness  matrix,  Kd  is  the  differential  stiffness  matrix,  and  the  h  are  the 
eigenvalues  to  be  computed.  Once  the  eigenvalues  are  determined,  the  critical  buckling 
load  is  found  from: 

P  =AP  (9) 

where  Pcr  are  the  critical  buckling  loads  and  Pa  is  the  applied  load. 

NASTRAN  incorporates  three  methods  of  eigenvalue  decomposition  for  buckling 
analysis  because  no  single  method  is  perfect  for  all  matrices.  As  the  size  of  the  matrix 
increases,  or  if  the  matrix  is  ill-conditioned,  the  manipulation  of  all  of  the  terms  in  the 
matrix  to  find  a  solution  leads  to  computer  round-off  errors  and  division  by  zero.  The 
methods  used  are  Lanczos,  Inverse  Power,  and  Sturm  Modified  Inverse  Power  [Ref.  5]. 


D.  CONVERGENCE  OF  ANALYSIS  RESULTS 

As  described  by  Bathe  [Ref.  6],  FEA  requires  the  idealization  of  actual  physical 
problems  into  mechanical  descriptions,  followed  by  finite  element  solutions  of  the 
idealization.  He  states,  “. .  .a  proper  FE  solution  should  converge  (as  the  number  of 
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elements  is  increased)  to  the  analytical  (exact)  solution  of  the  differential  equations  that 
govern  the  response  of  the  mechanical  idealization.”  Bathe  goes  on  to  highlight  the 
sources  of  error  in  a  FEA  which  include  “discretization,  numerical  integration  in  space, 
evaluation  of  constitutive  relations,  solution  of  dynamic  equilibrium  equations,  solution 
of  FE  equations  by  iteration,  and  round-off  error.”  [Ref.  6,  p.  165,  Fig.  4.23,  Table  4.4] 

In  light  of  this,  it  was  assumed  that  if  an  accurate  and  complex  mechanical  idealization  of 
a  real  strut  was  used  as  an  input,  and  a  FEA  converged,  an  accurate  solution  of  the 
buckling  load  should  result.  The  accuracy  of  the  mechanical  idealization  was  validated 
by  comparison  of  mass  properties.  It  was  hoped  that  the  sources  of  error  described  above 
would  be  covered  by  the  maturity  of  the  FEA  tool  in  use,  namely  NASTRAN. 

E.  MIDSIDE  NODES  AND  CONVERGENCE 

Two  types  of  nodes  used  in  the  definition  of  finite  elements  are  comer  nodes  and 
midside  nodes.  Comer  nodes  are  points  at  which  the  elements  connect  with  each  other. 
Midside  nodes  can  be  included  to  provide  interpolation  points  between  the  comer  nodes. 
The  convergence  on  a  solution  is  greatly  affected  by  the  presence  of  midside  nodes.  Two 
methods  of  convergence  used  in  FEA  are  polynomial,  or  P-convergence,  and  stepwise,  or 
H-convergence.  P-convergence,  through  the  use  of  midside  nodes,  includes  displacement 
interpolation  between  the  comer  nodes.  The  number  of  nodes  added  between  the  corners 
equals  the  order  of  the  interpolation  polynomial  plus  one.  Thus,  the  term  parabolic 
interpolation  is  used  when  one  midside  node  is  included  on  each  edge  of  the  element. 

This  is  a  result  of  the  works  of  B.  Szabo,  E.  Wilson,  and  K.  Bathe.  Stepwise,  or  H- 
convergence,  is  by  decreasing  step  size,  i.e.,  element  size.  It  is  computationally  efficient 
and  appropriate  for  large,  uncomplicated  structures. 

P-convergence  is  faster  with  respect  to  element  dimensions  and  computational 
speed  and  less  demanding  on  computational  resources  but  requires  more  complex 
mathematical  implementation  (Rayleigh-Ritz  within  each  element).  H-convergence  is 
slower  and  easier  to  implement,  but  much  more  computing  resource  demanding.  The 
current  trend  is  to  use  H-convergence  because  of  accessibility  of  increased  computational 
power.  Convergence  for  eigenvalue  decomposition  is  more  demanding  than  merely 
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computing  stresses  and  displacements  because  of  the  need  to  evaluate  all  possible 
displacements.  More  compliant  elements  converge  better  in  eigenvalue  problems 
because  of  the  smooth  transition  from  one  element  to  the  next.  However,  in  some  cases, 
compliant  elements  can  lead  to  erroneous  convergence  which  results  is  extremely 
convoluted  displacements.  An  example  of  this  can  be  seen  in  Figure  HI-1,  which  is  from 
one  of  the  aluminum  tube  buckling  models  discussed  in  Appendix  D. 


Figure  m-1.  Convoluted  Displacement  Results  From  NASTRAN  Buckling  Analysis 

F.  THREE  DIMENSIONAL  ELEMENTS 

The  type  of  element  used  in  this  analysis  was  a  solid  tetrahedron  with  four 
triangular  surfaces,  six  edges,  and  four  corner  nodes.  This  type  of  element  does  not  allow 
for  rotation  of  the  comer  nodes.  The  result  is  four  nodes,  with  three  translational  degrees 
of  freedom  each,  for  a  total  of  12  degrees  of  freedom  per  element.  Without  nodal 
rotations,  the  elements  tend  to  interlock  and  are  not  very  compliant  for  a  given  number  of 
degrees  of  freedom.  When  midside  nodes  are  included  on  the  six  edges,  the  total  degrees 
of  freedom  for  the  individual  element  increases  to  30.  Even  without  nodal  rotation,  the 
effect  of  the  midside  nodes  is  to  release  the  element  interlocks,  making  the  structure  more 
compliant.  The  two  different  kinds  of  solid  tetrahedron  elements  are  shown  in  Figure 
m-2  below. 
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Figure  HI-2.  Solid  Tetrahedron  Elements  With  &  Without  Midside  Nodes. 


G.  ASSUMPTION  OF  CONSISTENT  CONVERGENCE 

The  hypothesis  was  made  that  regardless  of  whether  or  not  midside  nodes  were 
present,  the  eigenvalue  decompositions,  if  they  converged,  would  converge  to  consistent 
solutions  between  models  with  the  same  element  characteristics.  This  allowed 
comparison  of  buckling  loads  among  different  configurations  of  drag  strut  corrosion 
models  with  the  same  element  type  and  size.  In  this  investigation,  the  corroded  model 
buckling  load  is  normalized  by  the  nominal,  non-corroded  model  buckling  load. 
Thereby,  any  convergence  difficulties  would  be  self-compensating.  Thus  the  percent 
reduction  in  buckling  load  from  corrosion  could  be  used  even  though  the  absolute  value 
of  the  load  may  not  be  precise. 
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IV.  MODELING 


A.  MODELING  BY  GEOMETRY 

The  modeling  process  began  with  the  original  Lockheed  technical  drawing  used 
to  manufacture  the  drag  strut  [Ref.  5].  The  computer-aided  design  (CAD)  software 
application  AutoCAD2,  Release  13,  was  used  to  create  a  three-dimensional  (3D)  solid 
representation  of  the  stmt  that  matched  the  engineering  drawing  nominal  dimensions.  All 
aspects  of  the  drawing  where  included  in  the  AutoCAD  model,  including  the  various 
chamfer  radii  at  surface  intersections.  The  details  of  the  modeling  process  are  included  in 
Appendix  C.  Figure  IV- 1  indicates  the  overall  flow  of  the  modeling  process  as  well  as 
the  high  level  of  replication  achieved.  The  goal  was  to  make  the  model  match  the  real 
stmt  as  closely  as  possible.  Several  protrusions  on  the  lug  ends  were  not  included  in  the 
model  because  they  would  have  only  increased  the  complexity  of  the  model  without  any 
contribution  to  the  strength  of  hollow  cylinder  portion  of  interest  in  this  investigation. 
However,  the  eccentricity  that  could  exist  in  a  manufactured  stmt  as  a  result  of  allowable 
tolerances  in  the  flash  welding  process  was  included.  Close  inspection  of  the  figure 
shows  the  presence  of  the  eccentricity  at  the  flash  weld,  indicated  by  the  arrow .  The 
AutoCAD  solid  model  was  exported  into  an  ACIS3  file,  which  is  one  of  many  CAD 
industry  standardized  file  formats  for  transporting  3D  solid  and  2D  surface  and  curve 
geometries  between  different  software  applications. 

B.  MODELING  BY  FINITE  ELEMENTS 

The  FE  modeling  process  involved  the  use  of  the  software  application 
MSC/NASTRAN  for  Windows4,  Version  3.0.2.  The  ACIS  file  from  AutoCAD  was 
imported  to  NASTRAN  as  a  single  3D  solid  geometry  along  with  all  the  curves  and 

2  AutoCAD  is  a  registered  trademark  of  Autodesk,  Inc. 

3  ACIS  is  a  registered  trademark  of  Spatial  Technology,  Inc. 

4  Windows  is  a  registered  trademark  of  the  Microsoft  Corporation. 
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AutoCAD  Model 


Figure  IV-1.  Modeling  Process  Flow  &  Sample  Drag  Strut. 


surfaces  that  defined  that  geometry.  The  built-in  NASTRAN  solid  meshing  routine  was 
used  to  define  a  grid  of  3D  solid  tetrahedron  FEs  with  four  comer  nodes  each.  The 
meshing  routine  could  also  create  tetrahedron  elements  using  six  midside  nodes  in 
addition  to  the  four  comer  nodes,  for  a  total  of  ten  nodes  per  element.  Both  types  of 
tetrahedron  elements  are  shown  in  Figure  ID-1.  The  impact  of  using  midside  nodes  in  the 
FEA  process  was  discussed  in  section  IE  above.  Midside  nodes  were  not  used  in  this 
investigation  because  the  license  purchased  for  MSC/NASTRAN  for  Windows  only 
permitted  models  of  up  to  5,000  nodes.  It  was  not  possible  to  model  the  stmt  with 
tetrahedron  elements  including  midside  nodes  with  only  5,000  total  nodes.  The 
maximum  element  size  parameter  was  varied  within  the  built-in  meshing  routine  until  a 
model  of  the  stmt  was  obtained  with  as  close  to  5,000  nodes  as  possible,  leaving  some 
room  for  model  growth  when  simulated  corrosion  effects  where  included. 
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C.  MODEL  VALIDATION  THROUGH  MASS  PROPERTIES 


Both  the  AutoCAD  and  NASTRAN  models  were  validated  in  terms  of  accurately 
representing  the  real  drag  strut  by  comparing  the  computed  mass  properties  to  measured 
values  from  the  sample  rejected  strut.  The  mass  properties  evaluated  were  the  total  mass, 
the  longitudinal  axis  center  of  gravity,  and  the  moment  of  inertia  about  one  of  the 
transverse  axes.  The  results  are  summarized  in  Table  IV-1  below  which  shows  that  both 
of  the  computer  models  match  the  sample  strut  to  within  five  percent.  In  regard  to  the 
mass  measurement,  the  lower  values  for  the  models  were  expected  because  the  synthetic 
rubber  coating  and  lug  protrusions  were  not  included.  The  detailed  procedure  used  to 
determine  the  measured  properties  is  included  in  Appendix  A. 


Property 

Measured 

AutoCAD 

NASTRAN 

Computed 

Percent 

Difference 

Computed 

Percent 

Difference 

Total  Mass  (lbm) 

28.61 

26.97 

-5.73 

27.21 

-4.89 

Center  of  Gravity  (in) 

13.38 

13.49 

0.82 

13.46 

0.60 

Moment  of  Inertia 
(lbm-in2) 

4027 

3935 

-2.28 

3986 

-1.02 

Table  IV-1.  Drag  Strut  Mass  Properties  Comparison 


D.  MODELING  OF  EQUIVALENT  SPHERICAL  PITS 

The  results  of  the  corrosion  attribute  measurement  described  in  Appendix  A  were 
used  to  create  a  model  that  included  all  of  the  grid  element  pitting.  A  cut-away  view  of 
half  of  the  cylindrical  portion  of  that  model  is  shown  in  Figure  IV-2  below.  The  pit 
modeling  was  performed  in  AutoCAD  as  described  in  Appendix  C.  The  pitting 
configuration  could  not  be  modeled  with  finite  elements  with  less  than  5,000  nodes  due 
to  the  level  of  geometric  complexity. 
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eluding  Grid  Element  Spherical  Pitting. 


30 


V.  DRAG  STRUT  BUCKLING  ANALYSIS  RESULTS 


A  total  of  34  drag  strut  configurations  were  modeled  and  analyzed  using  both 
AutoCAD  and  MSC/NASTRAN  for  Windows.  These  configurations  ranged  from  the 
nominal  strut,  with  no  corrosion  effects,  to  struts  with  several  two  inch  holes  through 
them.  The  analyses  included  the  linear  static  compression  case  and  two  to  three  linear 
buckling  cases.  In  addition,  several  linear  static  tension  cases  were  performed  to  uncover 
any  peculiarities  between  compression  and  tension.  As  discussed  above,  the  load 
boundary  conditions  in  tension  caused  an  incorrect  stress  concentration  in  the  lugs; 
however  the  far  field  stress  was  correct.  The  conclusion  was  that  the  highest  stress 
element  and  value  matched  the  compression  case  and  thus  tension  analyses  for  all  the 
configurations  was  unnecessary.  The  following  sections  contain  a  brief  description  of 
each  configuration,  along  with  a  diagram  and  the  buckling  and  yield  load  data.  All  of  the 
load  data  is  in  thousands  of  pounds  or  kilo-pounds  (kips).  The  nominal  case  was  used  to 
normalize  the  remaining  load  data  to  highlight  the  impact  of  the  defects  more  than  the 
absolute  significance  of  the  buckling  load  itself.  Therefore  all  of  the  other  configuration 
buckling  load  data  is  presented  as  a  fraction  of  the  nominal. 

A.  NOMINAL  CASE 

The  solid  geometry  created  in  AutoCAD  was  imported  into  NASTRAN  and 
meshed  with  finite  elements  to  produce  the  nominal  strut  without  corrosion  effects.  Two 
models  were  produced  from  this  configuration  with  different  maximum  element  sizes. 
The  first  had  an  element  size  of  0.415  inches  which  resulted  in  4,990  total  nodes  and 
14,723  elements.  This  was  below  the  5,000  node  limit  but  left  no  room  for  growth. 
Another  model  was  built  with  an  element  size  of  0.420  inches,  resulting  in  4,829  nodes 
and  14,247  elements.  This  became  the  standard  element  size  for  all  future  models  to 
assure  some  consistency  in  results.  The  result  of  a  buckling  analysis  on  the  nominal 
model,  including  the  deformed  shape  after  buckling  and  colored  Von  Mises  stress 
contours,  is  shown  in  Figure  V-l. 
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Figure  V-l .  Nominal  Drag  Strut  MSC/NASTRAN  Elastic  Buckling  Analysis  Results. 
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A  model  number  was  assigned  for  configuration  identification  among  all  34 
configurations  analyzed.  The  model  number  is  included  with  each  subsequent  strut 
image  for  identification.  The  buckling  analysis  on  the  nominal  strut  also  provided  the 
locations  for  the  critical  nodes  for  the  first  three  buckling  mode  shapes.  These  locations 
were  then  used  to  locate  pits  and  other  defects  in  future  configurations.  The  first  three 
buckling  mode  shapes  and  the  locations  for  the  critical  nodes  are  indicated  by  the  arrows 
in  Figure  V-2  below. 

1.  First  Buckling  Mode 

Buckling  Load:  Pcr=  729.0 

Node  Location:  x  =  15.083 
y  =  0.000 

z  =  1.888 

2.  Second  Buckling  Mode 

Buckling  Load:  Pcr  =  760.6 

Node  Location:  x  =  16.008 
y  =  -1.888 
z  =  0.000 

3.  Third  Buckling  Mode 

Buckling  Load:  Pcr=  2,215 

Node  Location:  x  =  5.598 

y  =  0.000 

z  =  1.500 

Figure  V-2.  First  Three  Buckling  Modes  (Units:  Loads  -  kips..  Locations  -  in.) 

B.  SIX  MOST  SEVERE  PITS 

The  goal  of  this  investigation  was  to  perform  an  analysis  of  the  lower  tail  of  the 
probability  distribution  of  the  population  of  all  struts.  Analyzing  the  worst  case  would 


33 


provide  an  upper  bound  for  the  risk  of  today’s  strut.  The  limitations  discussed  above 
prevented  analysis  of  a  model  that  included  all  of  the  corrosion  discovered  in  the  sample 
strut.  However  a  limited  model  was  analyzed  with  the  six  most  severe  pits  from  the 
sample  strut.  Two  models  were  constructed  with  two  different  maximum  element  sizes. 
The  added  complexity  of  the  pits  forced  maximum  element  sizes  of  0.430  and  0.423 
respectively.  Figure  V-3  shows  the  position  of  the  pits  and  Table  V-l  shows  the  load 
data. 
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Figure  V-3.  Six  Most  Severe  Pits. 


Element 

Yield 

Normalized 

Model 

Size 

Load 

Buckling 

Number 

(in) 

(kips) 

Load 

3 

0.430 

298.1 

0.995 

4 

0.423 

295.5 

0.994 

Table  V-l.  Six  Most  Severe  Pits  Load  Data. 


C.  COMPLETELY  THIN  SECTION 

As  an  extension  of  the  six  most  severe  pits  model,  a  model  was  created  with  an 
entire  sector  thinned  by  the  depth  of  the  deepest  pit  discovered,  or  0.041  inches.  This 
configuration  approximated  the  instability  in  the  pitted  sample  strut,  because  it  results  in 
a  model  that  is  quite  unbalanced  in  terms  of  local  centroid  and  moment  of  inertia.  Since 
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column  buckling  is  directly  related  to  the  stability  of  the  structure,  this  model  was 
assumed  to  be  as  unstable  as  any  pitting  that  would  be  found  in  an  operational  strut. 
Figure  V-4  shows  the  orientation  of  the  thinned  section  with  respect  to  the  strut.  The 
compressive  yield  load  was  249.9  kips  while  the  normalized  buckling  load  was  0.964. 


Section  2 A  Thinned  by  0.041  inches 


Figure  V-4.  Completely  Thinned  Section. 


D.  MORE  SEVERE  SPHERICAL  PITS 

The  next  step  was  to  find  the  upper  bound  of  the  risk  for  a  strut  more  corroded 
than  the  sample  strut.  This  would  represent  the  strut  of  the  future  that  is  allowed  to 
continue  to  corrode.  In  this  case  the  pits  were  deep  enough  to  extend  to  the  limit  of  the 
wall  thickness  or  0.160  in.  In  addition,  these  very  severe  pits  were  located  at  the  critical 
locations  determined  from  the  nominal  strut.  There  were  seven  models  created:  three 
configurations  of  one  pit  at  a  single  node,  three  configurations  of  combinations  of  two 
pits  at  two  different  nodes,  and  finally,  one  configuration  of  pits  at  all  three  nodes.  The 
spherical  pit  radius  was  set  at  twice  the  pit  depth  so  that  the  lip  of  the  pit  would  not  be 
perpendicular  to  the  wall.  This  was  done  to  match  the  majority  of  the  pits  observed 
during  the  sample  strut  dissection  and  is  explained  in  more  detail  in  Appendix  B.  Figure 
V-5  shows  a  cutaway  drawing  of  the  strut  with  the  three  possible  pit  locations  depicted. 
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Table  V-2  shows  the  various  combinations  of  pits  and  the  associated  yield  load  and 
normalized  buckling  load. 


Figure  V-5.  Spherical  Pit  Combinations  at  Nodes  1,  2,  &  3. 


Model 

Number 

Nodes 
With  Pits 

Yield 

Load 

(kips) 

Normalized 
Buckling 
Load  _ 

6 

1 

230.1 

0.996 

7 

2 

235.3 

0.999 

8 

3 

272.1 

0.998 

9 

1  &2 

230.5 

0.997 

10 

1  &3 

230.0 

0.994 

11 

2  &  3 

236.3 

0.997 

12 

1,  2  &  3 

229.9 

0.994 

Table  V-2.  Spherical  Pit  Combinations  Load  Data 

E.  HOLES  AT  NODE  1 

The  results  of  the  more  severe  spherical  pit  analysis  showed  extremely  small 
effect  on  the  buckling  load.  The  next  step  was  to  move  to  more  severe  configurations 
that  might  begin  to  reveal  the  point  at  which  pitting  could  become  critical  to  the  buckling 
failure  mode.  The  spherical  pits  were  replaced  with  through-the-thickness  holes  of 
increasing  diameter,  from  0.5  to  2.0  inches.  This  portion  of  the  analysis  was  confined  to 
node  1  because  it  is  the  most  critical  location.  Combinations  included  holes  at  node  1 , 
holes  on  the  opposite  side  of  node  1,  and  holes  on  both  sides  at  node  1.  Figure  V-6 
shows  the  various  hole  configurations  and  Table  V-3  contains  the  load  data. 
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Model 

Number 

Hole 

Size 

(in) 

Sides 

With 

Holes 

Yield 

Load 

(kips) 

Normalized 

Buckling 

Load 

13 

Front 

176.7 

14 

Back 

182.7 

15 

0-5 

Both 

182.6 

0.989 

16 

Front 

139.2 

0.981 

17 

1.0 

Back 

155.1 

0.981 

18 

Both 

145.6 

0.961 

19 

1.5 

Front 

113.1 

0.953 

20 

1.5 

Back 

98.8 

0.954 

21 

1.5 

Both 

95.5 

0.911 

22 

2.0 

Front 

93.1 

23 

2.0 

Back 

73.4 

24 

2.0 

Both 

64.6 

0.822 

Table  V-3.  Node  1  Hole  Combinations  Load  Data 
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F. 


HOLES  AT  NODE  2 


To  uncover  any  peculiarities  with  respect  to  buckling  mode  shapes  with  holes 
explored  above,  a  similar  but  less  exhaustive  approach  was  applied  at  node  2.  The  four 
configurations  analyzed  included  a  single  hole  at  node  2  of  diameters  from  0.5  to  2.0 
inches.  As  it  turned  out,  a  2.0  inch  hole  at  node  2  had  a  lower  buckling  load  than  the 
same  hole  at  node  1 ;  however,  this  was  not  the  case  with  the  yield  load.  Another 
peculiarity  that  the  node  2  analysis  revealed  was  that  the  first  buckling  mode  shape 
changed  when  a  1.5  inch  or  larger  hole  was  present.  The  shape  matched  that  of  the 
second  mode  in  a’1  prior  analysis.  Figure  V-7  shows  the  configurations,  along  with  the 
resulting  mode  shapes.  The  load  data  is  presented  in  Table  V-4. 


Figure  V-7.  Holes  at  Node  2  With  Mode  Shapes. 


Hole 

Yield 

Normalized 

Model 

Size 

Load 

Buckling 

Number 

(in) 

(kips) 

Load 

25 

0.5 

167.4 

0.999 

26 

1.0 

153.9 

0.997 

27 

1.5 

110.7 

0.990 

28 

2.0 

96.3 

0.876 

Table  V-4.  Holes  at  Node  2  Load  Data 
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G.  COMBINATIONS  OF  HOLES  AT  NODES  1  &  2 


Given  the  mode  shape  change  observed  above,  additional  combinations  of  holes 
at  both  nodes  1  and  2  were  analyzed  to  study  the  interaction  between  the  two  nodes.  In 
these  two  cases  a  1.5  inch  hole  was  placed  at  node  2  along  with  0.5  and  1.0  inch  holes  at 
node  1,  respectively.  The  case  with  the  0.5  inch  hole  at  node  1  still  exhibited  the  mode 
shape  change  observed  above;  however  the  1.0  inch  hole  forced  the  first  mode  shape 
back  to  the  original  case.  This  demonstrates  the  dominance  of  holes  (or  critical  pits)  at 
node  1.  Figure  V-8  shows  the  two  cases  analyzed  in  addition  to  the  resulting  mode 
shapes.  The  load  data  is  presented  in  Table  V-4. 


Figure  V-8.  Hole  Combinations  at  Nodes  1  &  2  With  Mode  Shapes 


Hole 

Yield 

Normalized 

Model 

Size 

Load 

Buckling 

Number 

(in) 

(kips) 

Load 

29 

Node  1  -  0.5 
Node  2- 1.5 

0.989 

Node  1-1.6 
Node  2  - 1.5 

86.4 

0.978 

Table  V-5.  Hole  Combinations  at  Nodes  1  &  2  Load  Data 

H.  OTHER  HOLE  COMBINATIONS 

All  the  preceding  configurations  contained  pits  or  holes  only  at  critical  nodes.  To 
verify  that  those  cases  are  the  worst  case  and  to  study  the  effect  of  clustering  holes 
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around  one  location,  the  following  configurations  were  analyzed.  The  first  was  a  0.5  inch 
hole  three  inches  closer  to  the  end  of  the  strut  than  node  1  (x  =  12  inches).  The 
expectation  was  that  this  would  yield  a  higher  buckling  load  than  the  case  of  the  same 
hole  at  node  1.  The  results  yielded  a  normalized  buckling  load  of  0.997,  which  is  higher 
than  the  value  of  0.993  found  model  number  10  above.  The  other  cases  analyzed 
included  clustering  of  three  0.5  inch  holes  around  node  1 .  In  one  case  the  holes  were 
spaced  1.0  inches  apart,  along  the  x  axis.  In  the  other  two  cases,  the  holes  were  arranged 
20  and  40  degrees  apart,  respectively,  around  the  circumference  at  node  1.  All  four  of 
these  cases  are  shown  in  Figure  V-9.  The  load  data  is  presented  in  Table  V-6. 


Figure  V-9.  Other  Hole  Combinations. 


Model 

Number 

Hole 

Size 

(in) 

Hole 

Location 

Yield 

Load 

(kips) 

Normalized 

Buckling 

Load 

31 

0.5 

x  =  12  in. 

180.7 

0.997 

32 

0.5 

Node  1 

1  in.  Apart 

195.6 

0.984 

33 

0.5 

Node  1 

20°  Apart 

120.1 

0.985 

34 

0.5 

Node  1 

40°  Apart 

180.7 

0.989 

Table  V-6.  Other  Hole  Combinations  Load  Data. 
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I.  BUCKLING  ANALYSIS  SUMMARY 

The  result  of  all  the  finite  element  modeling  and  analysis  is  the  effect  of  corrosion 
on  the  buckling  strength  of  the  strut  in  terms  of  normalized  buckling  load.  The 
conclusion  is  that  even  severe  corrosion,  much  more  severe  than  would  ever  be  accepted 
in  operational  use,  has  minimal  impact  on  the  buckling  strength.  This  is  highlighted  in 
Table  V-7  which  is  a  listing  of  all  the  strut  configurations  modeled,  sorted  by  corrosion 
severity. 
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Table  V-7.  Summary  of  Drag  Strut  Finite  Element  Analysis  Results. 
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VL  RISK  ASSESSMENT  FROM  A  SINGLE  RANDOM  VARIABLE 


The  risk  assessment  for  the  drag  stmt  in  the  presence  of  corrosion  was  made  by 
dividing  the  problem  into  two  distinct  modes.  The  buckling  risk  assessment  was 
performed  by  comparing  the  probability  of  occurrence  of  corrosive  severity,  mapped  into 
buckling  load,  to  the  operational  service  load.  The  yield  and  fracture  assessment  should 
be  performed  by  direct  filtering  with  a  low  temperature  proof  test. 

A.  COLUMN  BUCKLING  RISK  ASSESSMENT 

Certain  significant  corrosion  configuration  buckling  loads  were  determined  from 
finite  element  analysis  as  previously  described.  Each  of  the  32  cases  became  data  points 
in  a  probability  distribution  in  terms  of  buckling  load  and  the  probability  of  occurrence 
for  the  pitting  required  to  create  the  configuration  analyzed.  All  of  the  corrosion 
configurations  that  could  exist  in  the  stmt  could  not  possibly  be  modeled;  therefore,  a 
complete  probability  distribution  could  not  be  assembled.  However,  an  estimate  of  a 
probability  distribution  function  (PDF)  could  be  fit  on  top  of  the  32  data  points  that  were 
modeled.  The  relationship  of  corrosion  severity  to  buckling  load  is  neither  linear, 
monotonic,  nor  homologous  because  of  the  importance  of  the  location  of  the  pitting 
relative  to  the  critical  nodes  in  the  stmt.  This  is  demonstrated  in  Figure  VI- 1,  which 
shows  the  32  corrosion  configuration  buckling  loads  and  their  respective  joint  probability 
of  occurrence.  The  assumption  was  made  that  the  real  stmt  PDF  would  be  in  the  form  of 
a  Weibull  distribution  because  a  severe  pit  in  a  critical  location  would  be  the  weak  link 
that  would  cause  the  stmt  to  buckle.  A  Weibull  distribution  curve  fit  was  performed  until 
a  PDF  resulted  that  reasonably  encompassed  all  of  the  analysis  data.  The  shape  and  scale 
parameters  for  that  PDF  are  shown  in  Figure  VI- 1. 

With  the  shape  and  scale  parameters  from  the  PDF,  a  cumulative  distribution 
function  (CDF)  could  also  be  plotted  as  a  function  of  buckling  load.  The  CDF  provided 
the  crucial  relationship  necessary  to  perform  the  risk  assessment;  however  one  more  step 
was  necessary  to  make  the  CDF  relate  to  the  operation  envelope  of  the  drag  stmt.  The 
buckling  loads  were  determined  using  a  linear-elastic  buckling  analysis.  As  expected,  the 
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values  of  the  computed  loads  were  significantly  higher  than  the  design  load  or  the 
operational  load.  This  indicated  that  plastic  deformation  of  the  strut  must  occur  prior  to 
buckling.  The  discussion  included  in  Appendix  D  shows  the  significance  of  this 
phenomenon  demonstrated  through  buckling  testing  on  aluminum  tubes  that  matched  the 
buckling  characteristics  of  the  drag  strut.  The  relative  impact  of  corrosion  severity  on  the 
linear-elastic  buckling  solutions  was  assumed  to  be  transformable  to  the  operational  and 
design  load  regime  of  the  real  strut.  This  transformation  was  performed  by  adjusting  the 
scale  parameter  of  the  Weibull  distribution  to  that  of  the  manufacturers  estimate  of  the 
critical  buckling  load  for  the  nominal  drag  strut  [Ref.  1].  This  transformation  is  shown 
graphically  in  Figure  VI-2.  Included  on  the  figure  is  the  operational  envelope  of  the  strut 
in  compression  load,  displayed  at  the  Man-Safe  reliability  threshold  of  one  failure  out  of 
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Buckling  Load  vs.  Cumulative  Probability  of 
Occurrence 


_ Weibull  CDF  (Shape=100,  Scale=800) - Shifted  Weibull  CDF  (Shape=100,  Scale=278.5) 
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Figure  VI-2.  Buckling  Load  CDF  Transformation. 

100,000  struts,  or  a  probability  of  occurrence  of  1x10 5.  This  is  the  risk  assessment  for 
compressive  buckling  of  the  stmt.  The  shifted  Weibull  curve  represents  the  buckling 
load  given  the  joint  probability  of  occurrence  of  severe  pitting  at  a  critical  location.  The 
curve  indicates  that  it  is  extremely  improbable  that  corrosion  will  cause  buckling  in  the 
operational  load  environment  of  the  drag  stmt. 

B.  YIELD  AND  FRACTURE  RISK  ASSESSMENT 

The  risk  assessment  for  yielding  or  fracture  from  fatigue  could  be  performed  in  a 
traditional  manner  by  conducting  statistically  significant  testing  on  actual  hardware.  This 
method  is  indirect  because  it  does  not  assure  the  reliability  of  any  particular  stmt,  only 
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the  general  reliability  of  the  population  of  struts.  In  the  presence  of  corrosion,  this 
method  is  inadequate  because  one  could  never  achieve  statistical  significance  in  the  face 
of  the  corrosion  variability.  A  direct  method  for  not  only  assessing  the  risk,  but  assuring 
the  reliability  of  every  strut  is  a  low  temperature  proof  test. 

The  rationale  for  performing  a  low  temperature  proof  test  comes  from  estimation 
of  the  time-temperature  relationship  of  first  order  kinetics.  Every  10°  C  decrease  in 
temperature  is  equivalent  to  a  reduction  in  time  by  one  half.  If,  for  example,  a  proof  test 
were  performed  at  a  particular  load  for  100  seconds  at  a  temperature  of -200°  C  it  would 
be  equivalent  to  performing  the  same  test  for  0. 1  milliseconds  at  0°  C.  This  can  be  seen 
from: 

W,  =  f  ■ '  2"  =  (100  )ip*  )=  0.0001  sec  (10) 

This  time  is  too  short  for  any  flaw  growth  to  occur;  therefore,  the  proof  test  does  not 
cause  damage  to  any  strut  that  passes  the  test.  By  the  same  application  of  the  first  order 
kinetics  approximation,  any  flaw  that  already  exists  in  the  strut  will  cause  catastrophic 
failure.  The  time  shift  causes  the  loading  at  low  temperature  to  be  equivalent  to  an 
impact  load  at  normal  temperature  and  the  strut  behaves  as  if  it  were  extremely  brittle. 
Therefore  there  is  no  time  for  dislocations  to  move  and  no  yielding  occurs,  but  a  critical 
flaw  causes  failure.  The  result  is  a  direct,  high  pass  filtering  of  any  weak  struts,  with  no 
damage  to  the  surviving  struts.  The  test  could  be  performed  using  liquid  nitrogen  which 
has  a  boiling  point  of -196°  C  (-320°  F)  or  approximately  216°  C  below  room 
temperature.  Cooling  the  strut  to  liquid  nitrogen  temperatures  would  provide  the 
necessary  time  shift  as  shown  above.  The  low  temperature  proof  test  would  assure  that 
each  tested  strut  would  neither  yield  nor  fracture  in  fatigue  at  the  operational  load. 

It  would  not  be  necessary  to  perform  both  a  tension  and  a  compression  proof  test. 
Because  of  the  equivalence  of  the  tensile  and  compressive  yield  criterion,  such  as 
described  by  Von  Mises  stress  or  maximum  shear  stress,  it  is  only  necessary  to  perform  a 
tensile  proof  test.  If  a  strut  does  not  yield  in  tension,  then  it  would  not  yield  in 
compression  either. 
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To  provide  an  additional  level  of  confidence  in  the  fatigue  life  of  a  corroded  strut, 
several  struts  that  have  been  rejected  by  current  inspection  methods  could  be  proof  tested 
at  low  temperature.  The  ones  that  pass  could  then  be  fatigue  tested,  either  as  a  whole 
strut,  or  as  dog  bone  coupons.  The  data  from  these  tests  could  be  used  to  determine  the 
both  the  remaining  fatigue  life  and  the  residual  strength  of  the  surviving  struts.  This 
information  would  also  indicate  the  required  interval  for  subsequent  proof  tests  for  each 
specific  strut  after  continued  service. 

C.  COMBINED  RISK  ASSESSMENT 

The  two  individual  risk  assessments  combine  to  provide  the  quantitative 
assessment  necessary  to  assure  corroded  strut  reliability.  The  buckling  analysis  showed 
the  following  key  points: 

1.  Severe  corrosion  at  a  critical  location  is  highly  improbable.  The  joint 
probability  of  occurrence  of  a  severe  corrosive  pit  at  a  specific  location  is  less  than  one  in 
a  million  or  lxlCT6. 

2.  Column  buckling  of  the  strut  occurs  only  after  plastic  yielding  has  occurred. 

3.  The  finite  element  solutions  provided  very  accurate  plastic  yield  load  data  and 
the  plastic  yield  loads  were  always  less  than  the  buckling  load  for  the  same  configuration. 

4.  A  corroded  strut  would  have  to  yield  prior  to  buckling. 

5.  The  transformed  CDF  indicates  that  buckling  due  to  corrosion  is  extremely 
improbable. 

If  a  particular  strut  survives  a  low  temperature  proof  test,  then  that  strut  will  not 
yield  at  the  operational  load  and  therefore  it  will  not  buckle.  In  addition,  the  strut  would 
not  possess  a  flaw  large  enough  to  be  critical  in  terms  of  flaw  growth  and  fatigue  fracture. 
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VII.  CONCLUSIONS 


An  investigation  was  performed  on  the  effect  of  internal  corrosion  on  the  buckling 
strength  of  P-3C  main  landing  gear  lower  drag  stmts.  Based  upon  the  results  of  both  the 
corrosion  measurement  and  probability  distribution  definition  process  and  the  finite 
element  analysis,  the  probability  of  a  drag  stmt  buckling  in  operational  service  was  found 
to  be  less  the  one  chance  in  a  million,  which  is  extremely  remote.  Therefore  the  column 
buckling  failure  mode  for  a  severely  corroded  P-3C  main  landing  gear  lower  drag  stmt  is 
not  an  issue. 

A  low  temperature  proof  test  could  be  performed  that  would  assure  that  a  specific 
stmt  would  neither  yield  nor  fracture  at  the  operational  load.  The  tensile  yield  test  would 
also  assure  that  compressive  yielding  would  not  occur  because  of  equivalent  yield 
criterion.  Additionally,  if  the  stmt  would  not  yield,  then  it  could  not  buckle  either. 

The  reliability  assurance  for  the  drag  stmt  in  the  presence  of  corrosion  is 
summarized  in  Figure  VII- 1.  The  operational  environment  is  contained  within  the 
filtered  domain  defined  by  the  low  temperature  proof  test.  The  hash-marked  portion  of 
the  compression  operating  region  represents  the  loading  condition  that  can  only  result 
from  full  braking  while  taxiing  in  reverse.  Because  the  fleet  operators  never  perform  this 
operation,  the  likelihood  of  the  stmt  buckling  is  even  more  improbable.  The  finite 
element  analysis  of  buckling  load  region  is  separated  from  the  filtered  domain  by  the 
assurance  domain.  The  size  of  the  assurance  domain,  although  displayed  in  load,  is  in 
terms  of  probability  of  occurrence  of  either  corrosion  severity  or  the  time  required  to 
produce  the  corrosion. 
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Reliability  Assurance  Summary 
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Figure  VII- 1.  Reliability  Assurance  Summary. 
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VIII.  RECOMMENDATIONS 


Based  upon  the  results  of  this  investigation,  the  following  recommendations  are 
provided  for  the  P-3  Class  Desk  Officer,  Program  Manager  Air,  Maritime  Surveillance 
Aircraft  (PMA-290): 

1.  Perform  a  low  temperature  proof  test  on  main  landing  gear  lower  drag  strut 
currently  in  service.  Place  the  surviving  struts  back  in  service  for  a  length  of  time 
commensurate  with  the  confidence  level  of  the  fatigue  life  remaining.  Determination  of 
that  fatigue  life  was  beyond  the  scope  of  this  investigation. 

2.  Develop  a  method  to  prevent  further  corrosion  by  either  sealing  the  internal 
cavity  of  the  strut  and  filling  it  with  an  inert  gas,  such  as  nitrogen,  or  coating  the  internal 
surface  with  a  corrosion  inhibitor.  Other  possibilities  include  placing  some  kind  of 
organic  or  inorganic  desiccant  in  the  cavity  to  remove  the  moisture  that  provides  the 
corrosive  environment.  The  analysis  in  this  investigation  suggests  that  continued 
corrosion  may  still  be  acceptable  in  terms  of  buckling  because  of  the  extremely 
improbable  effect  of  corrosion  on  buckling  load.  However,  the  effect  of  continued 
corrosion  upon  the  fatigue  life  of  the  strut  was  beyond  the  scope  of  this  investigation. 

3.  Perform  a  low  temperature  proof  test  on  several  previously  rejected  drag 
stmts.  Use  the  surviving  stmts  for  both  fatigue  life  testing,  either  at  the  coupon  level  or 
as  a  whole  stmt,  and  then  residual  strength  testing  to  determine  the  service  life  remaining 
as  discussed  above. 
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APPENDIX  A.  REJECTED  DRAG  STRUT  PROCESSING 


1.  MASS  PROPERTIES 

a.  Total  Mass 

The  mass  of  the  sample  drag  strut  was  determined  by  direct  measurement  on  a 
digital  scale  with  automatic  zeroing  capability.  The  mass  was  determined  to  be  28.61  ± 
0.01  pounds. 

b.  Center  of  Gravity 

The  center  of  gravity  along  the  longitudinal  (x)  axis  was  measured  by  first 
marking  the  strut  using  the  computed  value  from  the  AutoCAD  solid  model.  Then  the 
strut  was  suspended  with  fine  steel  wire  around  the  cylindrical  portion  at  that  location,  as 
shown  in  Figure  A-l.  The  position  of  the  wire  was  then  adjusted  in  the  x-direction  until 
the  strut  hung  with  the  longitudinal  axis  in  a  level  attitude.  The  new  wire  position  was 
marked  and  then  measured  relative  to  original  computed  position.  The  x-axis  center  of 
gravity  was  measured  at  13.38  inches  from  the  origin,  which  is  the  center  of  the  small 
hole  directly  next  to  the  large  lug  end. 


Figure  A-l.  Center  of  Gravity  Measurement. 
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c.  Moment  of  Inertia 

The  moment  of  inertia  about  the  transverse  (z)  axis  was  measured  by  setting  up  a 
compound  pendulum  problem  with  the  stmt  suspended  from  a  single,  nearly  frictionless, 
pivot  point.  The  stmt  was  supported  by  stainless  steel  wire  tied  around  the  lug  ends  as 
shown  in  Figure  A-2.  The  stmt  was  hung  so  that  the  longitudinal  axis  was  level  with  the 
pendulum  arrangement  at  rest.  The  center  of  gravity  position  was  verified  once  again  by 
observing  that  it  was  directly  below  the  pivot  point.  The  stmt  was  given  a  series  of  gentle 
pushes  until  it  was  observed  to  swing  in  a  steady  arc  with  rotation  purely  about  the  z  axis. 
The  period  of  oscillation  was  determined  by  noting  the  time  required  for  the  stmt  to 
complete  30  complete  oscillations.  This  timing  process  was  accomplished  a  twice  to 
verify  accuracy.  The  entire  experiment  was  repeated  a  second  time  with  the  stmt 
suspended  at  a  different  pendulum  length. 


Once  the  period  of  oscillation  was  found  for  the  particular  setup,  the  moment  of 
inertia  about  the  center  of  gravity  was  calculated  from  the  following  relationship: 


where  lcg  is  the  stmt  moment  of  inertia  about  its  center  of  gravity,  w  is  the  weight 
of  the  stmt,  /  is  the  length  from  the  pivot  to  the  center  of  gravity,  T  is  the  period  of 
oscillation,  and  g  is  the  acceleration  due  to  gravity.  This  formula  is  derived  from  the 
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relationship  between  the  period  of  oscillation  and  the  pendulum  length  for  a  simple 
pendulum  and  the  parallel  axis  theorem  for  transformation  of  moments  of  inertia  about 
one  location  to  another  location.  The  mass  and  moment  of  inertia  of  the  wire  were 
neglected  due  to  their  insignificant  contribution  compared  to  the  mass  of  the  strut. 

d.  Calculations 

The  following  MATLAB  script  file  was  used  to  calculate  the  moment  of  inertia 

from  the  collected  data: 

%  moment_i .  m 

%  Calculate  the  moment  of  inertia  of  an  object  about  its  own 
%  c.g.  given  the  period  of  a  compound  pendulum  test. 

format  compact;  clear  all? 

r  =  [12. 64; 12. 64; 33. 06; 33. 06] ;  %  distance  from  pivot  to  c.g.  (in) 

w  =  28.61;  %  wieght  of  the  body  (lbs) 

t  =  [46. 95; 46. 86; 58. 64; 58. 49] ;  %  time  for  30  oscillations  (sec) 

g  =  32.174;  %  acceleration  of  gravity  (ft/secA2) 

T  =  t/30 ;  %  Period  of  oscillation  (sec) 

leg  =  w*(r/12)  .*  (T.  A2/ (4*piA2)  -  (r/12)  /g)  ;  %  (slug-ftA2) 

Icg_lbm  =  Icg*g*144;  %  (lbm-inA2) 

disp ( '  ' ) 

disp(' Period  of  oscillation  in  seconds') 

disp(T) 

disp ( '  ') 

disp ('Moment  of  inertia  in  slug-ftA2') 
disp (leg) 
disp('  ') 

disp ('Moment  of  inertia  in  lbm-inA2') 
disp ( Icg_lbm) 

The  following  text  is  the  output  from  the  moment  of  inertia  script  file: 

%  OUTPUT  % 

%»  moment_i 
% 

%Period  of  oscillation  in  seconds 
%  1.5650 

%  1.5620 

%  1.9547 

%  1.9497 

% 

%Moment  of  inertia  in  slug-ftA2 
%  0.8830 

%  0.8758 

%  0.8785 
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%  0.8395 

% 

%Moment  of  inertia  in  lbm-inA2 
%  1.0e+003  * 

%  4.0910 

%  4.0579 

%  4.0700 

%  3.8894 

%» 

2.  DISSECTION,  CORROSION  REMOVAL,  AND  MARKING 

Once  the  mass  properties  were  determined,  the  strut  was  dissected  into  sections  to 
allow  inspection,  corrosion  removal,  and  marking  of  the  grid  elements.  The  strut  was 
first  cut  across  the  diameter  into  five  pieces:  the  two  tapered  lug  ends,  the  two  cylindrical 
mid-sections,  and  the  center  area  containing  the  flash  weld.  The  two  cylindrical  mid¬ 
sections  were  then  cut  in  half,  lengthwise,  to  allow  access  to  the  internal  wall  surface. 

The  resulting  four  semi-cylindrical  pieces  contained  the  sections  of  interest  for  the 
buckling  investigation  because  the  material  in  these  sections  carry  the  bending  stress 
when  the  strut  is  put  in  compression.  The  remainder  of  this  discussion  refers  to  these 
four  sections.  Figure  A-3  shows  the  results  of  the  dissection. 


Figure  A-3.  Drag  Strut  Following  Dissection. 


There  were  large  layers  of  scale  that  easily  flaked  off  the  surface  upon  initial 
inspection.  Some  of  this  scale  is  shown  in  Figure  A-4  below.  These  layers  were  most 
likely  a  combination  of  excess  material  from  the  flash  weld  process  and  the  synthetic 
rubber  coating  that  had  been  hardened  by  the  baking  process.  The  remainder  of  the 
corrosion  was  removed  by  a  combination  of  mechanical  scraping  and  glass  bead  blasting. 
It  did  not  appear  that  the  blasting  process  affected  the  base  metal  because  the  fine 
machining  marks  from  initial  manufacture  were  still  visible.  Finally,  the  sections  were 
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scribed  with  grid  markings  at  one  half  inch  by  20  degree  intervals.  Figure  A-5  shows  an 
example  of  a  clean  semi-cylindrical  section  with  the  grid  elements  scribed  in  the  surface. 


Figure  A-4.  Semi-Cylindrical  Section  With  Flaking  Scale. 


Figure  A-5.  Semi-Cylindrical  Section  With  Scribed  Grid. 
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3.  CORROSION  EFFECTS  QUANTIFICATION 

Once  the  corrosion  was  removed  from  the  strut  and  grid  established,  the  effects  of 
the  corrosion  process  were  quantified  as  follows: 

a.  Pit  Depth 

A  dial  gage  capable  of  measuring  linear  displacement  to  0.001  inches  was 
mounted  in  a  fashion  that  allowed  direct  measurement  of  the  section  wall  thickness  and 
pit  depth.  The  dial  gage  was  fitted  with  a  sharp  point  that  allowed  penetration  into  the 
narrowest  pits.  The  depth  of  the  deepest  pit  in  each  grid  element  was  measured  and 
recorded  as  shown  in  Table  A-l.  The  deepest  pit  was  0.041  inches  with  a  remaining  wall 
thickness  of  0. 1 19  inches.  Thirty-four  of  the  pits  were  deep  enough  that  the  remaining 
wall  thickness  was  less  than  the  rejection  criteria  of  0.140  inches.  T  wall  thickness 
without  corrosion  across  all  four  sections  varied  from  0.165  to  0.154  inches.  Figure  A-6 
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Table  A-l.  Grid  Element  Pit  Depth  (x  10'3  in.). 
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shows  the  grid  element  pit  depth  as  a  3D  surface  plot.  The  orientation  of  the  plot  is  90° 
to  the  right  of  the  table.  Figure  A-7  shows  the  pit  depth  measuring  set  up. 
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Figure  A-6.  Strut  Corrosion  Pit  Depth  Surface  Plot. 


b.  Pit  Volume 

The  next  characteristic  of  interest  was  the  volume  of  material  removed  during  the 
corrosion  process.  Several  measurement  schemes  were  considered  including  pressing 
clay  into  the  pits  and  then  measuring  the  weight  of  the  clay  after  removing  the  pieces 
from  each  grid  element.  Then  the  weight  could  be  converted  to  volume  with  the  known 
density  of  the  clay.  However,  the  method  of  choice  was  more  direct  and  involved  filling 
the  pits  in  a  grid  element  with  a  liquid  while  precisely  measuring  the  volume  of  liquid 
required  to  fill  the  pits.  Alcohol  was  chosen  as  the  liquid  because  it  spread  evenly  into 
the  pits  and  formed  a  level  surface  across  the  top  without  beading  up.  A  micro-pipette,  as 
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Figure  A-7.  Pit  Depth  and  Pit  Volume  Measuring  Equipment. 


shown  in  Figure  A-7,  was  used  to  dispense  the  alcohol  into  the  pits.  This  device  allowed 
volume  measurement  in  increments  as  small  as  0.001  milliliters  (ml)  which  corresponds 
to  6.102xl0'5  cubic  inches  (in3).  The  largest  volume  of  pits  in  one  grid  element  was 
measured  to  be  0.039  ml.  This  equates  to  approximately  300  pieces  of  0.5  millimeter 
(mm)  pencil  lead  or  about  the  size  of  a  kernel  of  com.  Grid  elements  with  pitting  smaller 
than  0.001  ml  were  considered  to  have  no  pitting  at  all.  The  pitting  volume  in  each  of  the 
grid  elements  is  presented  in  Table  A-2  below. 
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Table  A-2.  Grid  Element  Pit  Volume  (x  10‘3  ml), 
c.  Equivalent  Spherical  Pit  Radius 

Because  the  exact  shape  of  the  pits  was  not  critical  to  a  column  buckling  analysis 
there  was  no  attempt  made  to  represent  the  shape  attribute.  Instead  the  depth  of  the  pit 
and  the  volume  of  material  missing  were  chosen  as  values  to  be  represented.  A  method 
was  needed  to  allow  proper  inclusion  of  these  parameters  in  AutoCAD  and  NASTRAN 
3D  solid  geometries.  A  spherical  pit  was  chosen  as  a  good  compromise  between 
modeling  ease  and  the  desire  to  accurately  represent  as  much  corrosion  effect  as  possible. 
The  spherical  pit  size  needed  to  be  chosen  so  that  both  the  pit  depth  and  volume  were 
represented  simultaneously.  The  equivalent  spherical  pit  radius  was  defined  as  the  radius 
of  a  sphere  such  that  a  slice  from  that  sphere,  at  a  depth  equal  to  the  grid  pit  depth,  would 
occupy  the  same  volume  as  the  grid  pit  volume.  This  is  more  clearly  illustrated  in  Figure 
A-8  and  is  calculated  from  equation  A1  as  follows: 


63 


3V  I 

/?=-  +  (Al) 

3  ;r/i2 

V  y 

where  R  is  the  equivalent  spherical  pit  radius,  h  is  the  pit  depth,  and  V  is  the  pit  volume. 
A  justification  for  the  choice  of  this  attribute  in  terms  of  corrosion  pit  propagation  can  be 
found  in  Appendix  B.  The  equivalent  pit  spherical  radius  values  for  the  grid  elements  is 
presented  in  Table  A-3. 

Spherical  Radius,  R 


Pit  Depth,  hN 


,Pit  Volume,  V 


Cylindrical ^ZZZZ///ZZZZA 

Figure  A-8.  Equivalent  Spherical  Pit  Radius. 
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Table  A-3.  Grid  Element  Equivalent  Spherical  Pit  Radius  (in). 
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d.  Pit  Severity 

The  final  corrosion  attribute  used  in  this  investigation  was  the  pit  severity  which 
has  been  described  in  the  Probabilistic  Reliability  chapter  of  the  main  body.  It  was  not 
measured  directly  but  derived  from  the  combination  of  pit  depth  and  pit  volume.  Pit 
severity  was  chosen  for  the  development  of  a  probability  distribution  because  as  a  single 
attribute  it  combined  two  quantifiable  attributes  that  directly  effect  the  buckling  strength 
of  the  drag  strut.  The  computed  grid  pit  severity  values  are  shown  in  Table  A-4  below. 
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APPENDIX  B.  EQUIVALENT  SPHERICAL  PIT  JUSTIFICATION 


4.  GENERAL  CORROSION  PROCESS 

The  high  strength  steel  used  in  the  drag  strut  is  very  susceptible  to  corrosion  as 
evidenced  by  its  position  next  to  aluminum  and  magnesium  on  the  Galvanic  Series  for 
Seawater,  which  is  a  widely  used  ranking  or  approximation  of  how  metals  behave  in  a 
corrosive  environment.  Because  the  metal  is  in  a  highly  refined  form,  it  is  inherently 
unstable  and  is  continually  trying  to  revert  to  its  natural  mineral  state.  Unlike  some  other 
naturally  self-protecting  materials  such  as  copper,  corroding  steel  forms  a  loose,  porous, 
and  non-adhering  oxide  that  tends  to  flake  off  and  expose  fresh  metal  to  the  corrosive 
environment.  [Ref  8]  This  flaking  process  within  the  drag  strut  is  shown  in  Figure  A-4. 

5.  NUCLEATION  SITES 

The  corrosive  environment  is  established  when  moist  air  is  forced  into  the  strut 
during  the  numerous  decompression  and  compression  cycles  of  flight  operations.  It 
appears  from  observation  of  the  sample  strut  that  the  moisture  collects  or  beads  up  into 
droplets  in  certain  areas,  either  on  the  top  or  bottom  of  the  cylindrical  portion.  This  is  a 
naturally  occurring  phenomenon  relating  to  surface  energy  that  can  be  observed  on  the 
hood  of  a  car  after  a  rain.  These  droplets  form  localized  corrosive  cells,  or  nucleation 
sites,  in  which  the  corrosion  initially  propagates  at  the  same  rate  in  all  directions,  forming 
small  spherical  pits.  The  propagation  rate;  however  is  influenced  by  the  molecular 
characteristics  at  any  one  site,  which  accounts  for  the  rough,  irregular  boundaries  of  the 
corrosion  sites.  When  several  nucleation  sites  are  adjacent,  the  pits  eventually  grow 
together,  causing  sharp  boundaries  between  pits.  These  boundaries  have  higher  potential 
and  therefore  corrode  at  a  higher  rate.  Because  the  boundaries  between  pits  corrode 
faster  than  the  outside  edges,  the  many  small  spherical  pits  grow  together  into  one  big  pit. 
Also,  because  of  the  higher  propagation  rate  at  the  sharp  boundaries,  the  slope  at  the  edge 
of  the  pits  becomes  more  shallow  than  at  the  initial  nucleation  sites.  Figure  B-lshows  an 
idealization  of  the  process  described  above,  which  is  the  justification  for  the  equivalent 
spherical  pits  used  in  modeling  as  describe  in  Appendix  A. 
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Initial  nucleation  Omni-directional  Propagation 

(Auto-corrosion) 


Sharp  Boundaries  Equivalent  Spherical  Pit 


Figure  B-l.  Corrosion  Propagation  Idealization. 
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APPENDIX  C.  DRAG  STRUT  MODELING  PROCEDURE 


1.  AutoCAD  PROCEDURES. 

a.  Engineering  Drawing. 

The  engineering  drawing  used  to  manufacture  the  strut  was  closely  examined  and 
annotated  with  all  of  the  information  that  would  be  necessary  input  for  constructing  a 
CAD  drawing.  This  included  the  center  coordinates  for  all  chamfer  radii  or  other 
machining  processes.  All  coordinates  were  referenced  to  the  manufacturer’s  origin  which 
was  the  center  of  the  0.562  inch  diameter  hole  on  the  large  lug  end  of  the  strut.  The  x- 
axis  was  selected  to  run  down  the  longitudinal  centerline.  The  y-axis  was  selected  so  that 
the  centerline  of  the  0.562  inch  diameter  hole  would  be  perpendicular  to  the  x-y  plane 
with  the  large  lug  end  of  the  strut  on  the  left.  The  coordinate  system  has  been 
superimposed  on  the  engineering  drawing  as  shown  in  Figure  C-l. 

b.  X-Y  Plane  Cross  Section. 

The  cross  section  made  by  cutting  the  strut  through  the  center  with  the  x-y 
coordinate  plane  was  drawn  in  AutoCAD  using  a  series  of  line  segments,  arcs  and  circles. 
Only  the  upper  portion  of  the  cross  section  was  drawn  as  if  the  strut  were  cut  by  the  x-z 
coordinate  plane.  These  lines  were  then  trimmed  as  appropriate  to  form  smooth 
intersections  as  indicated  in  the  engineering  drawing.  The  result  was  one  closed  line  that 
described  the  exterior  profile,  the  interior  profile  of  the  cylindrical  portion,  and  the 
maximum  radius  of  the  large  and  small  lug  ends.  The  individual  line  segments  that  made 
up  this  closed  line  was  then  joined  together  into  one  continuous  polyline.  This  polyline 
was  then  revolved  360  degrees  about  the  x-axis  to  form  a  solid  cylinder  with  hollow 
center  and  lug  ends. 

c.  Virtual  Machining  and  Trimming  Process. 

The  machining  process  was  simulated  by  constructing  solid  primitives  (blocks, 
cylinders,  and  spheres)  and  locating  them  relative  to  the  solid  cylinder  from  above  so  that 
when  subtracted  from  the  drawing  the  required  solid  geometry  would  remain.  In  some 
cases  the  coordinate  system  had  to  be  tilted  first  in  order  to  properly  align  the  primitives 
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prior  to  subtraction.  The  result  of  this  process  was  a  virtual  strut  geometry  that  matched 
the  real  strut  as  described  in  Appendix  A.  Some  of  the  protrusions  on  the  lug  ends  of  the 
strut  were  not  modeled  as  they  would  only  serve  to  increase  the  complexity  of  the  model 
without  contributing  to  the  analysis.  It  was  expected  that  the  mass  properties  of  the 
model  would  then  be  slightly  less  in  magnitude  than  the  real  strut.  The  mass  properties 
data  computed  directly  by  AutoCAD  is  shown  below. 


P-3C  LOWER  DRAG  STRUT  AutoCAD  MODEL  ( strut_01 . dwg)  05/21/98 

Mass  Property  Computations  with  origin  at  0.562"  diameter  hole 
-  SOLIDS  - 


Mass : 

95.0449 

Volume: 

95.0449 

Bounding  box: 

X: 

-5.5750 

33.6800 

Y: 

-2.6350  -- 

2.6350 

Z: 

-2.6350  — 

2.6350 

Centroid: 

X: 

13.4858 

Y: 

-0.0013 

Z: 

0.0000 

Moments  of  inertia: 

X: 

210.3237 

Y: 

31099.6935 

Z: 

31153.4474 

Products  of  inertia: 

XY: 

-0.2087 

YZ: 

0.0024 

ZX: 

-0.0021 

Radii  of  gyration: 

X: 

1.4876 

Y: 

18.0890 

Z: 

18.1046 

Principal  moments  and  X-Y-Z  directions  about  centroid: 

I:  210.3234  along  [1.0000  0.0001  0.0000] 

J:  13814.0552  along  [-0.0001  1.0000  0.0000] 
K:  13867.8088  along  [0.0000  0.0000  1.0000] 


Mass  Property  Computations  with  origin  at  Centroid 
-  SOLIDS  - 


Mass : 

95.1037 

Volume : 

95.1037 

Bounding  box: 

X: 

-19.0533  — 

20.2017 

Y: 

-2.6337  — 

2.6363 

Z: 

-2.6350  — 

2.6350 

Centroid: 

X: 

0.0000 

Y: 

0.0000 

Z: 

0.0000 

Moments  of  inertia: 

X: 

210.3651 

Y: 

13822.8298 
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z 

13876.5692 

Products  of  inertia:  XY 

1.4484 

YZ 

-0.0007 

ZX 

0.0260 

Radii  of  gyration:  X 

1.4873 

Y 

12.0559 

Z 

12.0793 

Principal  moments  and  X-Y-Z  directions  about  centroid: 

I:  210.3649  along  [1.0000  0.0001  0.0000] 

J:  13822.8299  along  [-0.0001  1.0000  0.0000] 
K:  13876.5692  along  [0.0000  0.0000  1.0000] 

The  mass  and  moments  of  inertia  presented  above  are  in  terms  of  volume  and 
have  units  of  in3  and  in5  respectively.  They  had  to  be  converted  to  the  values  presented  in 
Appendix  A  by  applying  the  material  density  for  AISI 4340  steel  which  is  0.2836  lbm  per 
in3. 

d.  Eccentricity  Representation 

The  manufactured  drag  strut  can  have  a  certain  amount  of  eccentricity  as  a  result 
of  the  flash  weld  process.  The  engineering  drawing  allows  for  as  much  as  0.045  inches 
mismatch  in  the  diameters  of  the  two  cylindrical  sections  following  welding  and  heat 
treatment.  The  AutoCAD  model  was  sliced  by  a  y-z  plane  at  the  center  of  the  flash  weld 
and  the  right  half  was  displaced  0.0318  inches  in  both  the  positive  y  and  z  directions. 

The  two  halves  were  then  joined  together  again  into  one  3D  solid.  This  resulted  in  a  total 
of  0.045  inches  of  offset  causing  an  internal  moment  in  the  strut  when  loaded  axially. 

e.  Corrosion  Pitting  Representation. 

The  pitting  measured  during  the  procedures  outlined  in  Appendix  A  was  included 
in  the  AutoCAD  model  as  follows.  The  process  involved  creating  spherical  primitives  at 
specific  locations  within  the  strut  and  then  subtracting  those  spheres  which  left  a 
spherical  pit.  The  spheres  were  located  according  to  the  data  presented  in  Table  C-l 
which  is  the  pit  center  radial  distance  in  each  of  the  grid  elements.  The  left  most  column 
in  the  table  shows  the  longitudinal  position  and  the  top  row  shows  the  angular  position 
around  the  cylinder  for  each  grid  element  center.  First,  the  AutoCAD  work  plane  was 
rotated  about  the  x-axis  according  to  the  angular  position  for  a  column  of  data.  Then,  a 
sphere  was  created,  centered  at  the  longitudinal  position  in  the  x  direction  and  the  radial 
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distance  in  the  y  direction.  The  radius  of  the  sphere  was  set  according  to  the  data  in 
Table  A-3  which  is  the  equivalent  spherical  pit  radius  for  each  grid.  Then  this  sphere  was 
subtracted  from  the  strut,  leaving  the  spherical  pit.  An  example  of  the  pitting  is  shown  in 
Figure  C-2. 
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Figure  C-2.  Strut  Model  Including  Grid  Element  Spherical  Pitting. 
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f.  ACIS  File  Export 

Following  the  modeling  procedures,  the  AutoCAD  3D  solid  geometry  was 
exported  into  an  ACIS  file  to  facilitate  use  by  NASTRAN. 

2.  NASTRAN  FOR  WINDOWS  PROCEDURES 

The  following  is  a  list  of  the  abbreviated  steps  accomplished  in  NASTRAN  for 
Windows  for  each  of  the  34  configurations  of  drag  struts  that  were  modeled. 

a.  File  /  Import  /  Geometry . . . 

1)  Change  to  directory  where  ACIS  file  resides 

2)  Solid  Model  Read  Options:  OK  (Accept  default) 

3)  NASTRAN  loads  solid  geometry  with  all  curves  and  surfaces  that 
define  the  geometry. 

b.  Model  /  Material. . . 

1)  ID:  1 

2)  Load:  AISI 4340  Steel 

3)  Mass  Density:  change  to  0.2836  (So  that  mass  properties  come  out 
in  correct  units.) 

4)  OK /Cancel 

c.  Model /Property... 

1)  ID:  1 

2)  Title:  Drag  Strut  Solid 

3)  Elem/Property  Type:  Volume  Elements  /  Solid  /  OK 

4)  Material:  AISI  4340  Steel  (From  list) 

5)  OK  /  Cancel 

d.  Geometry  /  Curve  /  From  Surface  /  Slice. . . 

1)  Solid  ID:  1 

2)  Methods:  Coordinate  System  Plane  -  yz  plane  (x  =  -2.94) 

3)  Repeat  for  x  =  3 1 .97 

(This  step  is  necessary  to  add  two  curves  to  the  solid  geometry  that  will 
allow  future  selection  of  nodes  on  a  surface  for  applying  boundary  conditions.) 

e.  File  /  Save  As. . .  (Save  the  geometry  alone  using  a  descriptive  filename) 

f.  Mesh  /  Geometry  /  Solid . . . 

1)  Element  Size:  0.420  (Or  as  appropriate  for  model  complexity  to 
stay  under  5,000  node  limit) 

2)  OK 
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3)  Property:  l..Drag  Strut  Solid 

4)  Midside  Nodes:  OFF  (Unless  model  complexity  allows) 

5)  OK 

6)  If  5,000  node  limit  exceeded  then  Tools  /  Undo  (CTRL  +  z)  to 
remove  the  nodes  and  try  again  with  larger  maximum  element  size. 

g.  Modify  /  Renumber . . . 

1 )  Node  . . .  /  Select  All  /  X  coordinate  order 

2)  Element. . .  /Select  All  /  X  coordinate  order 

h.  Model  /  Node . . .  (Add  nodes  for  bearings) 

1)  Coordinates:  (-2.94, 0, 0) 

2)  Coordinates:  (31.97,0.0318,0.0318) 

i.  Model /Element... 

1)  Type:  Rigid  /  OK 

2)  Node  (Independent) 

3)  Select  bearing  center  node  ' 

4)  Nodes  (Dependent) 

5)  Method:  On  Surface 

6)  Select  bearing  surfaces  /  OK 

7)  OK 

8)  Repeat  for  other  bearing 

j.  Model  /  Load  /  Nodal. . . 

1)  Name:  Tension  or  Compression 

2)  Select  large  bearing  node 

3)  Force/  Vector  /Specify... 

4)  Change  to  snap  to  node  (Right  mouse  click) 

5)  Select  start  node 

6)  Select  finish  node 

7)  Magnitude:  1 

8)  OK 

k.  Model  /  Constraint  /  Nodal .. . 

1)  Name:  Bearing 

2)  Select  bearing  node 

3)  Select  DOF  (Large:  2/3/4,  Small:  1/2/3) 

4)  OK 

5)  Repeat  for  other  bearing 

l.  Customize  View 

1)  View  Quick  Options  (CTRL  +  q) 

2)  All  Entities  Off 

3)  Check:  Elements,  Loads,  Constraints 
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4)  Done 

5)  View  /  Select  (F5) 

6)  Full  Hidden  Line  -  OK 

m.  File  Rebuild  /  Yes  /  Yes 


n.  File  /  Save  As. . .  (Save  the  finite  elements  in  new  descriptive  name) 


o. 


P- 


File  /  Save  As... 

1)  Save  copy  in  F:\Scratch\*.mod  (This  forces  all  scratch  files  created 
during  analysis  to  be  placed  in  the  Scratch  directory.) 

File  /  Analyze... 

1)  Type:  Buckling  (or  Static  for  Tension) 

2)  Output  Types:  Displacements  and  Stresses 

3)  Eigenvalues:  2  or  3  (Buckling  only) 

4)  OK 

5)  OK  to  save  file. . .  Yes 


3.  MAKING  PITS  WITH  NASTRAN 


a.  View  /  Options. . .  (F6) 

1)  Tools  and  View  Style 

2)  Workplane  and  Rulers 

3)  Check  Draw  Entity 

4)  OK  (Turns  the  workplane  and  rulers  on  so  one  can  verify  that  the 
workplane  is  properly  aligned  for  future  coordinate  references.) 

b.  View  /  Rotate. . .  (F8)  -  As  required  so  that  the  required  x-y  plane  will  be 
in  the  plane  of  the  screen  with  the  z-axis  coming  out  of  the  screen. 

c.  Tools  /  Workplane  /  Select  Plane. . . 

1)  Method... /Align  to  View 

d.  View  /  Rotate. . .  (F8)  -  As  required  to  allow  3D  verification  of  solid 
primitive  direction  and  placement.  Trimetric  works  well.) 

e.  Geometry  /  Solid  /  Primitives  /  Sphere. . . 

1)  Radius:  0.318  (or  as  desired) 

2)  Coordinates:  (Example  for  0.318  radius  sphere) 

Model:  15.083,-1.5690 

Mode  2:  16.008,-1.5372 
Mode  3:  5.598,-1.1580 

3)  Remove  (Negative  direction) 
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4)  OK 

(Nodes  for  modes  1  &  3  direct.  Node  2  adjust  y  for  offset.) 

4.  MAKING  HOLES  WITH  NASTRAN 

a.  Perform  steps  3  a  through  3  d  above. 

b.  Geometry  /  Solid  /  Activate . . . 

1)  Select  the  desired  solid  (l..Drag  Strut  Solid) 

c.  Geometry  /  Solid  /  Primitives . . . 

1)  Select  Cylinder 

2)  Material:  Remove 

3)  Direction:  As  required  to  orient  the  solid  properly  (Use  the 
reference  arrow  as  a  guide.) 

4)  Origin:  Enter  x  and  y  coordinates 

5)  Primitive:  Select  Cylinder 

6)  Radius:  as  desired 

7)  Height:  as  desired  (will  build  from  the  x-y  plane  in  the  z  direction 
that  is  selected  as  shown  with  the  arrow.) 

8)  OK 
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APPENDIX  D.  ALUMINUM  TUBE  BUCKLING  TESTS 


1 .  TUBE  BUCKLING  TESTS 


The  absolute  accuracy  of  the  FE  buckling  solutions  for  the  drag  struts  could  not 
be  validated  by  an  actual  buckling  test  with  a  real  strut;  however  buckling  tests  were 
performed  on  straight  hollow  aluminum  tubes  with  fixed-end  conditions.  A  NASTRAN 
model  was  created  that  matched  the  geometry,  material,  and  loading  boundary  conditions 
of  the  tubes.  Because  both  the  tube  and  strut  models  matched  their  real  counterparts,  the 
buckling  analysis  results  could  be  compared.  A  search  was  made  of  available  stock  to 
match  the  critical  buckling  characteristics  of  the  strut  as  closely  as  possible.  The  first 
characteristic  was  the  diameter  to  thickness  ratio,  which  for  the  drag  strut  is  23.59.  The 
closest  match  was  found  at  19.23  with  seamless  A1 2024-T3  tubes  that  were  1.25  inches 
in  diameter  with  a  wall  thickness  of  0.065  inches.  Next,  the  length  of  the  tubes  was 
chosen  to  match  the  slenderness  ratio  of  the  drag  strut  (effective  length  divided  by  cross 
section  radius  of  gyration).  This  value  was  27.28  for  the  drag  strut  which  resulted  in  a 
test  length  of  1 1.5  inches  for  the  tubes.  Additional  tube  lengths  of  8,  10,  12,  13,  and  23 
inches  were  tested  to  provide  further  resolution  on  the  FE  results.  Finite  element 
solutions  were  computed  for  all  of  the  test  lengths,  including  additional  lengths  of  30  and 
45  inches  to  investigate  buckling  behavior  in  the  purely  elastic  region  of  the  material. 

The  test  and  analysis  results  are  presented  in  Table  D-l  and  Figures  D-land  D-2  below. 
In  addition  to  the  test  and  FE  solutions,  the  results  from  the  appropriate  form  of  the 
classical  Euler  column  buckling  equation  for  fixed  end  conditions  are  provided.  The 
form  of  the  Euler  equation  is  [Ref.  9]: 


,  _47tEI 
E  _  L2 


(Dl) 


The  results  show  that  the  FEA  buckling  load  values  and  the  Euler  equation  values 
were  significantly  higher  than  the  tested  values,  as  expected;  however  the  compressive 
yield  values  matched  very  closely.  The  buckling  mode  shape,  half  sinusoidal  wave, 
matched  what  was  observed  which  lends  credence  to  the  boundary  conditions  included  in 
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Aluminum  Tube  Buckling  Test  Data  -  8/1 3/98  &  9/3/98 


Sample 


Buckling 
Load,  Pcr 
(lbs) 


14200 


13450 


13500 


13150 


13050 


13000 


12900 


12850 


12650 


12750 


12750 


12300 


12500 


10400 


Mean 

Euler 

FEA 

Per 

Load,  PE 

Load,  Pfe 

(lbs) 

(lbs) 

(lbs) 

14200 

281173 

118010 

13475 

179951 

114742 

12933 

136068 

108932 

12750 

124966 

101734 

12400 

106480 

89187 

10400 

34017 

32094 

19995 

19333 

8886 

8765 

Al  2024-T3  Modulus,  E  =  1.0700E+07 


Tube  Moment  of  Inertia,  I  =  0.0426 


Notes 


Added  dial  gage  to  detect  lateral  deflection 


Began  strain  control  process  in  final  stage  of  loading 


Pcr  determined  from  extrapolation 


Buckling  Load  Data 

Fixed  End  Boundary  Conditions 

Seamless  Al  2024-T3  Tube,  1.250"  OD  x  0.065"  WT 

o  Test  Data  A  FEA . Power  Fit  Extrapolation - Euler 
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Figure  D-2.  Test,  Analysis,  and  Euler  Buckling  Data  Combined. 


the  model.  The  FE  solution  indicated  that  the  first  yielding,  on  an  individual  element 
level,  would  occur  around  10,000  pounds.  The  8  inch  test  length  was  selected 
specifically  to  determine  the  yield  load.  The  specimen  was  short  enough  to  fit  into  a 
digital  caliper  with  a  resolution  of  ±  0.0005  inches.  The  length  prior  to  test  was  measured 
at  8.0550  inches.  The  tube  was  placed  under  compressive  loads  in  500  pound  increments 
with  length  measurements  between  each  step  until  a  permanent  change  in  length  was 
observed.  This  occurred  at  a  load  of  10,000  pounds  when  the  measured  length  changed 
to  8.0545  inches. 
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2. 


INTERPRETATION  OF  RESULTS 


The  most  important  result  of  the  tube  buckling  tests  was  observation  of  the 
convergence  upon  accurate  FEA  solutions  when  the  physical  geometry  caused  buckling 
within  the  linear  elastic  region  of  the  material.  This  should  be  self  evident  because  the 
FEA  solution  was  performed  under  the  linear-elastic  assumption.  However,  this  testing 
revealed  the  order  of  magnitude  of  divergence  and  mono  tonic  relationship  between  the 
solution  and  the  real  test  when  the  linear-elastic  assumption  was  violated.  The  results 
indicate  a  predictable  trend  in  solution  values  as  the  geometry  of  the  tube  was  changed. 
The  conclusion  was  that  the  same  assumption  could  be  made  about  the  drag  strut  as 
physical  geometry  was  changed  due  to  corrosive  pitting.  In  addition,  the  results  indicate 
that  high  confidence  could  be  placed  in  the  determination  that  buckling  of  the  drag  strut 
could  only  occur  following  plastic  deformation.  Therefore,  because  the  FEA  solutions 
provided  a  very  accurate  prediction  of  yield  load,  if  one  could  assure  that  a  corroded  strut 
would  not  yield  under  an  operational  load,  then  the  same  assurance  could  guarantee  that 
it  would  not  buckle  under  that  load. 
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